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ABSTRACT 

We present a study on the shells (and bubbles) in the Perseus molecular cloud using the COM- 
PLETE survey large-scale 12 CO(1-0) and 13 CO(1-0) maps. The twelve shells reported here are spread 
throughout most of the Perseus cloud and have circular or arc-like morphologies with a range in ra- 
dius of about 0.1 to 3 pc. Most of them have not been detected before most likely as maps of the 
region lacked the coverage and resolution needed to distinguish them. The majority of the shells are 
coincident with infrared nebulosity of similar shape and have a candidate powering source near the 
center. We suggest they are formed by the interaction of spherical or very wide-angle winds powered 
by young stars inside or near the Perseus molecular cloud — a cloud that is commonly considered 
to be mostly forming low-mass stars. Two of the twelve shells are powered by high-mass stars close 
to the cloud, while the others appear to be powered by low or intermediate mass stars in the cloud. 
We argue that winds with a mass loss rate of about 10 -8 to 10 -6 M yr _1 are required to produce 
the observed shells. Our estimates indicate that the energy input rate from these stellar winds is 
similar to the turbulence dissipation rate. We conclude that in Perseus the total energy input from 
both collimated protostellar outflows and powerful spherical winds from young stars is sufficient to 
maintain the turbulence in the molecular cloud. Large scale molecular line and IR continuum maps 
of a sample of clouds will help determine the frequency of this phenomenon in other star forming 
regions. 

Subject headings: ISM: bubbles — stars: winds, outflows — stars: formation — stars: premain- 
sequence — turbulence — ISM: individual (Perseus) 



1. INTRODUCTION 

Stars go through a mass loss stage during the early 
stages of their life. The properties of the wind driven 
by the forming star will depend, in part, on the mass 
and evolutionary stage of the driving source. In general, 
young embedded protostars drive powerful collimated 
(jet-like) outflows, and as protostars evolve the mass loss 
rate decreases and the outflow becomes less c ollimated 
( |Bontemps et al ||1996| |Arce fc Sargent ||2QQ6[ ). During 
the late stages ol star formation, low and high mass pre- 
main sequence stars (T-Tauri and Herbig Ae/Be stars) 
drive wide-angle winds, that in some cases coexist with a 
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collimated jet-like wind component. Regardless of their 
degree of collimation, these supersonic winds can entrain 
and accelerate the surrounding gas, thereby injecting mo- 
mentum and energy into the surrounding gaseous envi- 
ronment. 

Early on in the study of winds from young stars it was 
realized that they have the potential to significantly af- 
fect the dynamics and structure of their parent cloud. 



Norman & Silk (1980) proposed a model in which the 
collision between bubbles formed by the interaction of 
spherical winds from low-mass pre-main sequence stars 
and the surrounding quiescent gas would dominate the 
dynamics of molecular clouds. Although the idea that 
winds from young stars ca n have a significa nt impact on 



the cloud is still relevant (lArce et al. 



2007), the bubble 



model of Norman & Silk ( |1980p has lost its original ap 



peal as T Tauri stars are rarely observed to have a major 
impact on the surrounding molecular gas. Instead, it is 
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the bi-polar collimated winds from younger embedded 
protostars that have been commonly observed to have 
a significant effect on their surroundings through well- 
known molecular outflows. 

Parsec-scale circular cavities (typically referred to as 
bubbles) are regularly found in regions of high-mass star 
formation, and are lik ely created by spherical winds from 
high-mass stars (e.g.,| Church well et al.| 2QQ6, 2007 Beau 
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mass stars evolve much taster than low-mass stars and 
reach the main sequence stage while they are still accret- 
ing material and are embedded in their natal cloud. Dur- 
ing the main sequence phase, high-mass stars can drive 
spherical winds with mass loss rates of ~ 10 -6 M yr _1 
that can easily create the observed bubbles. The conven- 
tional thought has been that in regions where predomi- 
nantly low and intermediate mass stars form such bub- 
bles do not exist, as the spherical winds from pre-main 
sequence stars are thought to be too weak to produce 
them. In this study we hope to cast away this notion. 

Here we present a study of circular shells detected in 
CO line maps of the Perseus molecular cloud complex, a 
cloud generally thought to predominantly form low-mass 
stars. We find twelve shells in the molecular gas that 
are very likely produced by powerful spherical (or quasi- 
spherical) winds from young stars. Ten of these are from 
presumed low or intermediate mass stars (thought to be 
late B or of later spectral type) within the cloud and 
two are from high-mass stars close to the cloud. We find 
that these shells have the potential to have a significant 
impact on the cloud, as they have enough energy to drive 
the turbulence in Perseus and to disperse a small fraction 
of the cloud's gas. Yet, the few shells that we find are too 
sparse to drive the shell-shell collisions and produce the 
cl umps predicted by the original T Tauri bubble model 



of Norman fc Silk| ( jl980[ ). 

The Perseus molecular cloud complex is a chain of 
clouds with a total mass of about 7 x 10 3 M^, and en- 
compasses a tota l area of about 70 pc 2 (|Evans et al.|2009| 
|Arce et aL]|2010| ). There is a span of distance estimates 
to Perseus that range from 230 pc to 350 pc, in part 
due to different parts of the cloud being at slightly dif- 
ferent distances along the line of sight and also due to 
the uncertainties in distance measurements (see Enoch 
et al. 2006, and references therein for discussion). To 
simplify our calculations, we adopt a single fiducial dis- 
tance of 250±50 pc, similar to other recent cloud-wide 
studies of Perseus. The Perseus cloud complex contains 
two rich protostellar clusters, IC 348 and NGC 1333, and 
a number of other regions that are mostly forming low- 
mass stars, where most young stars associated with the 
cloud have spectral type of A or later. There are however, 
three known B stars that reside in the cloud: HD 281159, 
a B5 star in the IC 348 cluster; BD+30 549, a B8 star in 
NGC 1333; and SVS 3, a B5 star in NGC 1333. As we 
will show in this work (and others have also suggested) 
there are at least two nearby high-mass stars that inter- 
act with the cloud through their winds. 

The entir e Pers eus region was first surveyed in 12 CO 
by Sargent ( |1979[ ), and since then it has been mapped 
in CO at different angular resolutions (all with beams 
> 1.3') by a number of authors (e.g., Bachille r fc Cer- 
nicharo 1986; Ungere chts fc T haddeus 1987; Pad oan et 
al.| |1999| |Sun et al.||2006| ). The COmaps used in this 



study (see below) cover approximately the same area as 
previous studies of the Perseus cloud, but at a signifi- 
cantly higher angular resolution. It is possible that many 
of the shells discussed here were not previously detected 
due to the low angular resolution of previous maps. Re- 
cently, the Perseus molecular cloud complex was also ob- 
served in its entirety in the mid- and far-infrared as part 
of the "From Molecular Cores to Plane t-forming Disks" 
(a.k.a., c2d) Spitzer Legacy Project (| j0 rgensen et al. 
2006} |Rebull et al.|[2007| |Evans et aLpSM] ). These data 
have an angular resolution that is significantly less than 
1', and many of the shells detected in the molecular line 
maps are also observed as circular nebulosities in the c2d 
infrared images. 

2. DATA 

We use the 12 CO(1-0) and 13 CO(1-0) line maps of 
Perseus, collected as part of the COMPLETE (Coordi- 
nated Molecular Probe Line Extinction Thermal Emis- 
sion) Survey of Star Forming Region^] described in de- 
tail in |Ridge et al] ( |2006a| ). The 12 CO J=l-0 (115.271 
GHz) and the 16 CO J=l-0 (110.201 GHz) line maps were 
observed between 2002 and 2005 using the 14-meter Five 
College Radio Astronomy Observatory (FCRAO) tele- 
scope with the SEQUOIA 32-element focal plane array. 
The dual-IF digital correlator allowed the simultaneous 
observation of 12 CO and 13 CO, and a total bandwidth 
of 25 MHz with 1024 channels in each IF, providing a 
velocity resolution of approximately 0.064 km s _1 and 
0.066 km s _1 in the 12 CO and 13 CO maps, respectively. 
The telescope beam at these frequencies is about 46", 
and the resulting maps are Nyquist-sampled with 23" 
pixels. All figures shown here present the 12 CO and 13 CO 
in the T\ scale. To convert to the main beam antenna 
temperature scale, the spectra at 115 and 110 GHz need 
to be divided by 0.45 and 0.49, respectively. The median 
1 — a rms per channel in the 12 CO and 13 CO maps are 
approximately 0.23 and 0.11 K, respectively, in the Tj[ 
scale. 

3. SHELL AND SOURCE IDENTIFICATION 

In the molecular line maps, the shells in Perseus are 
seen as bright rings or arc structures surrounding re- 
gions of less (or no) emission and many are coincident 
with infrared nebulosity seen in the MIPS 24 /am c2d 
map. These shells, which we denote as COMPLETE 
Perseus Shell (CPS) 1 through 12, are listed in Table [TJ 
Their positions in the cloud are shown with respect to 
the 24 /im emission and the 13 CO (1-0) integrated inten- 
sity emission of the Perseus region in Figures jl] and [2] 
These figures serve as general finding charts that show 
the overall distribution of the shells in the Perseus cloud 
complex. More detailed figures that show the emission 
and velocity structure of each CPS are shown below. 

The shells in Perseus were found through different 
methods. Some were found searching for high-veolcity 
features with a 3D visualization technique (se e Figu re |3l), 
as par t of the outflow study conducted by |Arce e t "al. 
( |2010[ ). Others were found by comparing the molecular 
gas distribution with the c2d IRAC and MIPS images in 
regions close to bright, circular and arc-like, nebulosities. 
A few were detected by studying how the morphology of 



1 see http://www.cfa.harvard.edu/COMPLETE 
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the gas changed with velocity in our molecular line data 
cubes (see below for more detail). The sample of shells 
listed here is most likely incomplete, as we probably fail 
to detect small shells with radius of less than 1.9' (~ 0.14 
pc) since it is hard to identify a shell structure with a di- 
ameter that extends less than 5 pixels in our map. 

The regions associated with the CO shells were easily 
identified as circular (or arc) structures brighter than 
the surrounding molecular gas. We used the CO maps 
to obtain an estimate of the radius, thickness and center 
of each shell. We created CO integrated intensity maps 
over the velocity range in which the shell was clearly 
detected (shown in Figures |7| to [29]) . The center of each 
CPS was first estimated by visually fitting a circle to the 
CO integrated intensity map of the shell. Using circles at 
different positions and with different radii, we estimate 
the error in the center position to be about ±45" for the 
shells that have a structure that delineates more than 
half a circle's circumference (i.e., CPS 1, 2, 3, 6, 9, 10, 
11, 12). In shells that show only an arc in CO, but 
exhibit an almost complete circular structure in the IR 
maps (i.e., CPS 4, 5) we used the 24 fim map to estimate 
the center position of the CPS (which gives similar error 
as using the CO maps). The position and size for shells 
that show arcs that only trace less than half a circle's 
circumference are less constrained. For these shells (CPS 
7, 8) the errors in the center position estimate are larger 
by a factor of four. 

We obtained the shell radius and thickness by fitting 
a gaussian to the azimuthally averaged profile of the CO 
integrated intensity map of each CPS (as shown in Fig- 
ure Hi. The position of the profile's peak gives the shell 
radius and the full width at half maximum (FWHM) 
gives the thickness of the shell. We used the shell radius 
and thickness to define the area of the annulus used to 
estimate the mass, hereafter 



4.1 



for 
to 



"mass region" (see 
more detail). Although the same procedure was usee 
obtain the mass region for each CPS, for some shells we 
had to constrain it to a smaller area to avoid including 
emission from unrelated structures (i.e., clouds, filaments 
or other shells) in the mass estimates (see discussion for 
each CPS below). 

In order to exclude ring structures that might simply 
be due to random patterns in the cloud (e.g., due to 
turbulence or overlapping filaments), we only considered 
circular or arc-like structures that have a velocity distri- 
bution corresponding to that of an expanding bubble or 
thick ring. For an expanding bubble fully immersed in 
a cloud and detectable at all velocities, the radius of the 
ring seen on the plane of the sky will increase as the ve- 
locity increases from the bluest velocity where the shell is 
detected (i.e., the front cap of the bubble) up to the ve- 
locity of the bubble source (i.e., the bubble's central LSR 
velocity), and will then decrease with increasing veloc- 
ity up to the reddest velocity where the shell is detected 
(i.e., the back cap of the shell). A model of an expanding 
bubble inside a turbulent cloud with uniform densit^] is 



'' This model has 8 parameters: 77,, dr^, Poffset, ^offset^ dr c \ l 



eyn, urrn.s, t jl u 

|Cazzolato fc Pineaultj ([2005]) . The parameters , dr^ , and V e xp 
describe trie radius, thickness and expansion velocity of the bubble, 
respectively, and p ffset an d v ff set are the offset in the bubble's 
central position and velocity with respect to the cloud's center 
and velocity. The ambient cloud is modeled as a face-on slab of 



and 0. The exp anding bubble is modeled following 



shown in Figure |5j In Figure [5^i, the integrated inten- 
sity map (over a limited velocity range close to the cloud 
velocity) shows a clumpy circular structur^] The cor- 
responding position-velocity (p — v) diagram that cuts 
through the center of the bubble is shown in Figure ^jp. 
Here the bubble is identified as a circular or ring-like 
structure where the highest velocities (along the line of 
sight) are seen at the center of the bubble, and the gas 
with positive (negative) velocities correspond to the back 
(front) hemisphere of the bubble. If molecular gas does 
not surround the entire bubble, then the gas displaced by 
the bubble (i.e., the shell) will be observed over a more 
limited range of velocities. For example, if the center of 
the bubble is displaced towards the back of the cloud, 
the gas will only show the front hemisphere of the bub- 
ble. In this case the p — v will resemble Figure^, where 
the bubble appears as a "V" or "U" structure. Likewise, 
for a bubble with the center displaced toward the front of 
a thin cloud only the back hemisphere will be observed. 
This will be seen as an upside down "V" or "U" struc- 
ture in the p — v diagram, as in Figure [5]i. The shells 
in Perseus were not detected at all velocities (similar to 
those shown in Figures [5fc and [5]i) which indicates that 
most (if not all) of the bubbles (or winds) responsible for 
the observed shells are not fully immersed in the cloud 
(see discussion in § 4.2). We thus considered shells that 
exhibit a detectable change in the radius with velocity 
consistent with the expected kinematics of an expand- 
ing bubble that may or may not be fully surrounded by 
molecular gas. 

We searched for a candidate source for each COM- 
PLETE Perseus Shell using SIMBAD, and the 2MASS 
and Spitzer c2d catalogs. We used these catalogs to look 
for young (pre-main sequence) stars associated with the 
cloud, close to the center of each shell. We expect that 
young stars (with age ~ 1— 2 Myr or less) will be the most 
likely drivers of winds with enough ma ss lo ss rate to drive 
the observed shells (see discussion in § 4.3). The velocity 
of young (low-mass) stars with respect to their ambi- 
ent cloud is expe cted to be relatively low, between 0. 1 
and 0.2 km s" 1 ( |Hartmann||2002| |j0rgensen et al.||2007[ ). 
Thus, at the distance of Perseus, we expect sources that 
are about 1 Myr of age to be within 3' from the center 
of the circu lar cavity delinea t ed by the shell. Using the 
c2d survey, J0rg ensen et al.| ( |2006[ ) show that the den- 
sity of YSO candidates (YSOc) in Perseus outside the 
IC348 and NGC1333 clusters is 43.3 YSOc deg" 1 . If we 
assume that the c2d YSOc catalog lists 90% of the total 
PMS in the cloud (i.e., the catalog is ~ 10% incomplete, 
see Evans et al. 2009), the density should then be closer 
to 48 young stars per square degree. Hence, there is only 

thickness dr c i ou d with uniform density. All of the material evacu- 
ated from the bubble interior is redistributed uniformly throughout 
the shell. The velocity field of the ambient medium has a disper- 
sion of Vrms- Its spatial structure is fractal- like, with a power- 
spectrum of P(k) ~ k~P ( |Stutzki et al.|l998| >. These parameters 
are used to generate three-dimensional density and velocity fields, 
which are then transformed into position-position-velocity coordi- 
nates to simulate observations. In Figure[5|the integrated intensity 
greyscale is proportional to the gas column density. 

3 In this model, which uses a uniform density cloud, the clumpy 
structure of the shell is solely due to the cloud's turbulence. We 
would expect that a shell in a real turbulent molecular cloud with 
non-uniform density would be as (or more) clumpy than the model 
shell presented here. 
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a 30% probability of finding a pre-main sequence star by 
chance within a circular area with a radius of 3', for shells 
outside the central region of the IC 348 and NGC 1333 
cluster. For small shells inside these clusters (i.e., CPS 

3, 9, 10, 11), we search for possible sources within 1-1.5' 
radius (depending on the shell size). In these two clus- 
ters the surface density of YSOc is relatively large and 
a few to several tens of pre-main sequence stars may lie 
within a radius of 1 to 1.5'. However, for all these small 
shells we can identify the most likely candidate source 
from the list of stars in the region. We find a realistic 
candidate source for nine of the twelve shells, which we 
list in Table [ll Their properties are given in Table [3] 

The ring-like (or arc-like) morphology of the Perseus 
shells, their velocity-dependent structure, and in (most 
cases) the clear association with a source or an IR neb- 
ulosity with a similar morphology indicate that most, if 
not all, of these shells likely trace the region where (very 
wide angle or spherical) winds fr om y oung stars inter- 
act with the molecular gas (see § |4.3[ ). It is clear that 
not all shells show the same kind oi evidence supporting 
their presumed wind-driven nature. We quantify this by 
assigning a "confidence level" to each CPS based on the 
pieces of evidence that support the hypothesis that these 
structures are wind-driven, as opposed to, for example, 
being chance superposition of cloud filaments or random 
structures formed by the cloud's turbulence. The levels 
range from 2 to 5, where the higher the grade the more 
supporting evidence exist of that particular shell. The 
grade is assigned by adding the total number of pieces of 
evidence for each shell, as shown in Table [4] 1) change in 
shell radius with velocity; 2) coincidence with IR nebu- 
losity of similar shape and size; 3) CO or IR emission that 
clearly delineates more than half a circle; 4) & p — v dia- 
gram that is consistent with that of an expanding bubble; 
and 4) a realistic candidate source near the shell center. 
We note that even shells that do not show a clear p—v di- 
agram consistent with that of an expanding bubble show 
evidence of changing radius with velocity, but their p — v 
diagrams are contaminated by nearby shells (as is the 
case for CPS 7) or are just limited by the fact that the 
CO shell does not trace a full circle (as in the case of CPS 

4, 5, 7 and 8). We believe that all the shells discussed 
here are likely to be wind-driven structures, even shells 
with a confidence level of 2. A lower grade simply means 
that there is less evidence in the current data to support 
the shell authenticity compared to shells with a higher 
grade. 

Below we give a description for each CPS, the area and 
velocities associated with the CO shell, and discuss their 
candidate powering source. For each CPS we present 
two figures. In one we present channel maps of the re- 
gion surrounding the shell. In the other figure we show, 
on the left panel, an integrated intensity contour map 
of the shell (in most cases overlaid on the MIPS 24 /urn 
image of the region), and on the right panel we present 
a p — v diagram of the shell. We use information from 
the integrated intensity map and the p — v diagram to 
fit a simple spherical expan ding bubble model (based on 
Cazzolato fc Pineault|[2005 ) to the shell structure in the 
channel maps. The model consists of three parameters: 
bubble radius, expansion velocity, and central LSR ve- 
locity. The radius is set by our estimate of the CPS 
radius (shown in Table [T]). We use the p — v diagram of 



each source to obtain a first guess of the central veloc- 
ity (e.g., the upper part of the "U" structure in Figure 
[5fc) and the expansion velocity of the bubble (e.g., the 
aifference between the central velocity and the lowest 
point of the "U" structure in Figure [5fc). We then fit by 
eye the expanding bubble model to the channel maps by 
slightly changing our first guess to the expansion velocity 
and central velocity, while maintaining the bubble radius 
constant. The uncertainty in the estimate of the bubble 
central LSR velocity is approximately ±0.5 km s _1 . In 
most of the shells the expansion velocity is not well con- 
strained as the CO shell that traces the expanding bubble 
is observed over a limited range of velocities. Hence, the 
derived expansion velocity is a lower limit. 

3.1. CPS 1 

CPS 1 is located about 17' west of the central part of 
the NGC 1333 cluster and it is only distinguishable in the 
13 CO map (see Figured). It was first detected by looking 
at the molecular line cnannel maps. The shell structure 
is most clearly seen at 7.1 < Vlsr < 7.7 km s _1 , where 
it shows an approximately circular morphology with ~ 3' 
(0.2 pc) thick walls (see Figure [6]). Yet, in the channel 
maps and p — v plot emission that might be associated 
with the expanding shell can be seen over a wider range 
of velocities, from Vlsr ~ 5.5 km s _1 to ~ 8.0 km s _1 
(see Figure [6] and Figure [7]). The morphology of the 
13 CO emission within this range of velocities exhibits a 
structure consistent with an expanding bubble with only 
the near-side hemisphere observable; the shell radius in- 
creases with increasing Vlsr an d the p—v diagram shows 
a "U" or "V" shape (similar to the model shown in Fig- 
ure]^). From the p — v diagram and the channel maps 
we estimate an expansion velocity of ~ 2.5 km s _1 , and 
a central velocity at Vlsr — 8-0 km s _1 . To the north of 
CPS 1 there is a filament, seen in Figure [6] at Vlsr > 7.5 
km s _1 , which is not related to CPS 1. To the east 
of CPS 1 lies gas associated with the central region of 
the NGC 1333 cluster, most prominent at Vlsr — 7.9 
km s _1 . Unlike most of the other shells, CPS 1 is not 
easily discernible in 12 CO, as for this region (and ve- 
locities) the 12 CO (1-0) emission practically fills the ring 
structure seen in 13 CO. Neither the IRAC or MIPS im- 
ages show any evidence of a ring-like structure in this 
area. 

SIMBAD lists only two stars projected within a 3' ra- 
dius of the center of CPS 1. Out of these, a likely source 
is the source SSTc2d J032740. 5+311540 (also known as 
2MASS J03274053+31 15392). This source is classified 
as a star with an IR excess starting at 24 //m. This is 
very likely a young star, but it was not classified as a 
YSO candidate as this area of the sky was not covered 
by the c2d IRAC maps and IR AC detection is need ed for 
a positive YSOc classification ( jEyans et al.|j2QQ9| . It is 
very likely that this c2d source is BD +30 543, which has 
been classified as an F2V star and has been esti mated 
to be at the same distance (220 pc) as NGC 1333 ([Cer] 
nis| |1990). There is no bright star at the (low- accuracy) 
position quoted in SIMBAD for BD+30 543, and the 
closest bright c2d and 2MASS point source is SSTc2d 
J032740.5+311540. 

3.2. CPS 2 
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This shell is located just south of CPS 1 in the western 
edge of the NGC 1333 region. Similar to CPS 1, the 
circular structure is only clearly discernible in the 13 CO 
map; it is barely seen in the 12 CO map, and there is 
no IR nebulosity in the Spitzer images associated with 
this shell. CPS 2 was first recognized when we studied 
the distribution of the 13 CO (1-0) emission as a function 
of velocity. The roundish structure is most clearly seen 
in at Vlsr between 6.8 and 7.2 km s _1 , but arcs that 
appear to be part of the expanding shell can be seen at 
LSR velocities ranging from about 5.8 to 8.0 km s _1 (see 
Figure [8]). The channel maps and p — v diagram of this 
shell are consistent with that of an expanding bubble 
with only the front (blueshifted) hemisphere. Similar to 
CPS 1, the p — v diagram shows a "U" shape with a 
minimum LSR velocity of ~ 5.5 km s _1 and top velocity 
of ~ 8.0 km s _1 (Figure [9|. From the p — v diagram and 
the channel maps we estimate an expansion velocity of 
2.5 km s _1 and the central velocity to be at Vlsr — 8.0 
km s _1 . The mass region for this source was constrained 
to avoid including emission associated with CPS 1. 

It is not clear which star powers CPS 2. SIMBAD 
shows only one star within 4.8' of the presumed shell 
center: VSS IX-51. This star has a very low measured 
optical polarization of 0.78 ± 0.37, which cou ld suggest 
that this is a foreground star ( Vrba et al. 1976 ). This area 
of the sky was not covered by the c2d 1RAC maps and 
thus no c2d YSO candidate lies within the circular struc- 
ture. There are, however, two point sources with sig- 
nificant 24/im emission (SSTc2d J032739.1+310516 and 
J032742. 1+310529) and nineteen 2MASS sources with 
reliable photometry (from the 2MASS point source cat- 
alog) that lie within 2.5' of the center. None of these 
stars are relatively more likely to be a candidate shell- 
powering source; any of these could be the candidate 
driving source of CPS 2. 

3.3. CPS 3 

BD+30 549 is the B8 star illuminating the NGC 1333 
optical nebula. The shell associated with this source was 
discovered by examining the 12 CO channel maps, where 
we detect a bright structure that has the characteris- 
tics of a small expanding bubble. The 12 CO data show 
strong emission in a region of about 2' in size near the 
position of BD+30 549, at V LS r ~ 7.2 km s _1 , with 
a peak intensity that is more than three times the av- 
erage intensity of the surrounding area (see Figure 10). 
As the LSR velocity increases, the morphology of the 
high intensity emission changes to a ring-like structure 
that slightly increases in size, consistent with the veloc- 
ity structure expected for an expanding bubble with its 
front cap seen at Vlsr ~ 7 km s _1 . The ring structure 
is discernible up to Vlsr ~ 8.2 km s _1 , where the shell 
walls show a peak intensity that is about a factor of ap- 
proximately 2 higher than the surrounding emission. At 
larger Vlsr, CO cloud emission is present in this region, 
but shell walls are not discernible. Our simple expand- 
ing bubble model indicates that there should be a ring 
like structure for Vlsr > 8.4 km s _1 . This might not be 
detected clearly in our data due to our relatively large 
beam (compared to the size of the shell) and intervening 
(optically thick) CO emission. Higher angular resolution 
observation of an optically thin line tracer is needed to 
fully understand the kinematics of this bubble. The p — v 



diagram shows a clear "V" shape with a minimum and 
maximum LSR velocity at ~ 6.8 and 8.2 km s _1 , respec- 



tively (Figures 11). From the channel maps and p — v 
diagram we estimate an expansion velocity of 1.2 km s _1 
and the central velocity to be at Vlsr = 8.1 km s _1 . 

The main beam temperature (T m b) at the front cap 
(at Vlsr °° 7.2 km s _1 ) is about 45 K (the highest in 
the entire Persues 12 CO map), and in the ring structure 
(at Vlsr ~ 7.9 km s _1 ) the peak T m 5 is about 30 K. 
Assuming the 12 CO is optically thick we can use 

T = ^ (1) 

ex ln[l + 5.53/(T peafe +0.82)] V ; 

to obtain an estimate of the excitation temperature 
(T ex ), where T pea k is the peak intensity corrected for 
main beam efficiency (Rohlfs & Wilson 1996). We find 
that the excitation temperature in these regions is about 
50 and 35 K, similar to that found in the molecular gas 
of high- mass star forming regions. It is very likely that 
such high temperature is due to the strong radiation from 
BD+30 549. 

The candidate source of this expanding bubble lies 
about 80" to the west of the presumed center of the shell. 
The bubble's asymmetry, with respect to the position of 
BD+30 549, could be due to a slightly higher density 
in the gas immediately west of the bubble, where there 
is strong 13 CO emission. There is prominent nebular 
emission at optical and IR wavelengths associated with 
this source that fills the region enclosed by the CO shell 
(see Figure 11 ). 



3.4. CPS 4 

The presumed center of this shell lies about 27' east of 
the group of protostars in the Bl cloud (see Figure 
and it is clearly noticeable as a circular structure of about 
26' in diameter in the Spitz er IRAC 5.8 and 8.0 /im and 



MIPS 24 an d 70 /im maps ( |j0rgensen et al J2006I [itebull 
et al.|2007 ). The existence oi the molecular gas shell was 
first noticed when we conducted a search for high veloc- 
ity 12 CO using th e 3D rendering technique described in 



Arce et al.| (120101) , and found high- velocity gas close to 
the IR nebulosity. Unlike the c2d IR images, the molecu- 
lar line emission does not show a clear ring structure. In 
this region, the gas associated with CPS4 is only clearly 
discernible in the northern and eastern part of the struc- 
ture. In the west the shell structure is confused with gas 
associated with the Bl cloud. The southern part of CPS 
4 is not covered by the molecular line maps (see Fig- 
ures [l2] and 13). 12 CO emission presumably associated 
with CPS 4 is first seen in at Vlsr ~ ~ 1 km s _1 in a 
small area in the northern edge of the circular IR nebu- 
losity, and as the LSR velocity increases more emission 
is seen along both the northern and eastern parts of the 
nebulosity (see Figure 121. At Vlsr > 3.0 km s 1 gas as- 
sociated with the Bl cloud "contaminate" the region and 
it is hard to discern the emission associated with CPS 4. 
The p — v diagram shows a clear velocity gradient in the 
CO arc defined by CPS 4, with a minimum and maxi- 
mum LSR velocity at ~ — 1 and 4 km s _1 , respectively 
(Figures 13). The CO velocity structure of the arc is 



consistent with it being part of the backside (redshifted) 
hemisphere of an expanding bubble. Using the channel 
maps and p — v diagram we determine that a reasonable 
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fit to the data is an expanding bubble with an expansion 
velocity of 5.0 km s _1 and a central velocity at Vlsr — 
-1.1 km s _1 . 

A SIMBAD search shows HD 22179 as the only star 
within a 3.8' radius of the estimated center of CPS 4 
that could be a candidate driving source of this shell. 
HD 22197 is a rel atively young star, with a reported age 
of 30 to 100 Myr ( |Hillenbrand et al.l|2008| |Roccatagliata| 
et al. 2009), is a r elatively bright X -ray source ( IRAS 
J03353(O+311336, |Voges et al.|1999] ) ? and exhibits mid- 
IR excess emission associated with a debris disk — all 
properties that make it a suitable candidate source for 
CPS 4. Howeve r, the estimated distance to this source 



is 100 ± 20 pc ( [Hillenbrand et al . 
et al.||2009|), significantly smaller 



the Perseus dark cloud. Although different sources give 
slightly different spectral types for HD 22179, all are 
within a few sub-classes and all agree that HD 22197 is 
an early G-type star (GO by SIMBAD; G3V by Guillout 
et al. 2009; G5IV by Li & Hu 1998). Using the photome- 
try of this source (reported in SIMBAD) and its spectral 
type provides an estimate of the dista nce consistent with 
that given by Hillenbr and et al.] (120081 ). In addition, from 
its kinematics |Guillout et al = ~ j|2009| ) estimated that HD 
22197 has a 30% probability ofTJemg part of the IC 2391 
association and a 50% probability of being a member of 
the Pleiades, two clusters that are significantly closer to 
the Sun than the Perseus molecular cloud. It is there- 
fore very likely that HD 22197 is not associated with the 
Perseus cloud. 

We searched for other possible sources for CPS 4 in the 
c2d catalog, but it shows no YSO candidate within a ra- 
dius of 3' from the center of the shell. There are, however, 
more than 1154 point sources in this area listed in the 
c2d point source catalog. Near the center of the shell the 
IRAC1 image shows a group of stars clustered around 
a relatively bright one, SSTc2d J033525.4+310925 (at 
about 36" from the center of CPS4) classified as "star" 
with the slope of the SED between between 2 and 20 
fim (or a) equal to -2.55, and near to it (at approxi- 
mately V away from the center of CPS 4) lies SSTc2d 
J033526. 2+310902, a source that is classified as "ris- 
ing", with a = 0.77. These two seem the most likely 
candidates from the long list of c2d sources in this re- 
gion. However, there are no distances to these sources 
and no way to determine if any of these are associated 
with (or are close to) the Perseus cloud. Note that 
the fact that both the central velocity of the shell is 
Vlsr ~ — 1 km s _1 , significantly different from the Vlsr 
of the cloud in this region (~ 7 km s _1 ), and that CPS 4 
seems to be tracing only part of the backside hemisphere 
of an expanding bubble, indicates that the source of this 
shell lies foreground to the cloud. In summary, either 
CPS 4 is driven by a yet undetermined young star close 
to the front of the Perseus cloud or it is a shell completely 
unrelated (and foreground) to the Perseus cloud driven 
by HD 22197. 

3.5. CPS 5 

The infrared nebulosity and shell associated with CPS 
5 is clearly evident in all Spitzer IRAC and MIPS bands, 
as wel l as in IRAS images of the region (|Ridge et al. 



(~ 2.8 pc), has been previously reported by several stud- 
ies (see Bally et al. 2007, and references therein) and all 
agree that HD 278942, at the center of the ring, is the 
star powering CPS 5. There is some disagreement on 
the spectral type of HD 278942, but most estimates in- 
dicate that it is a late O or early B star (see Bally et 
al. 2007 fo r more discussion). The study by Ridge et 
al. (2006b) indicates that the shell is mostly behind the 
Perseus cloud, and that it is interacting with parts of the 
molecular cloud. In the molecular line maps, this inter- 
action is most clearl y se en in the western part of the shell 
(see Figures 14 and|15|). 

We first noticed the molecular gas shell when, search- 
ing for outflows in the region, we detected high veloc- 
han the distance to ity 12 CO in the southwestern edge of the infrared shell 



2008 Roccatagliata 



using t he 3D visualization technique described in [Arce 



2006bJ |j0rgensen eFluT1|2006l |Rebull et al.p007[ ). This 



et al. (2010). CO emission associated with CPS 5 is 
observed over a range of velocities, from Vlsr ~ 2 to 
7 km s -1 , but the arc-like CO structure is most clearly 
seen at Vlsr = 5.5 km s _1 , where it shows a morphol- 
ogy very similar to the IR nebulosity (Figures [l4| and 
15). The peak CO emission in this arc (T m 5 ~ 22 K) is 
significantly brighter than the surrounding emission, by 
about a factor of two. Assuming the emission is opti- 
cally thick, we estimate the CO excitation temperature 
in this region to be about 25 K (using Equation 1). The 
increase in excitation temperature in this region is more 
likely due to the interaction between the wind from HD 
278942 and the cloud. At Vlsr > 6 km s" 1 CO emission 
from CPS 7, CPS 6 and the L1468 cloud start appearing 
in the region, and at Vlsr > 6.5 km s _1 there is little or 



no CO emission associated with CPS 5 (see Figures 14). 

The p — v diagram for CPS 5 shows a clear velocity 
gradient in the CO emission. In an expanding bubble, 
this velocity structure is most likely to arise from an arc 
that is part of the front hemisphere (consistent with HD 
278942 being background to the cloud). The limited ex- 
tension of the CO emission (i.e., an arc that is much less 
than half the circumference of the ring traced by the IR 
emission) makes it hard to constrain the shell's expansion 
velocity. We use V exp = 6 km s _1 and a central velocity 
of Vlsr = 6.8 km s _1 for CPS 5 as it results in a rea- 
sonable fit of the expanding bubble model to the channel 
maps (Figure 14) and it is consistent with the limited 
p — v diagram oFthis source. We note that the expansion 
velocity of this source could easily be higher given that 
we do not to trace CO emission associated with this shell 
for angular offsets between and 30' in the p — v diagram 
shown in Figure 15 The mass region for this source is 



prominent circular structure, with a radius of about 38' 



constrained to the part of the CPS 5 annulus that lies in 
the right ascension range between 03^36 m and 03^40 m , 
and declination from 30°54 / to 32°08 / , in order to avoid 
emission from the L1468 cloud. 

3.6. CPS 6 

About 18' due east of HD 278942 lies IRAS 
03382+3145, an infrared source that is located at the cen- 
ter of a ring of about 5' in radius, seen in Spitzer MIPS 
maps ( jRebull et al. 2007). At near- to mid-IR wave- 
lengths, the MSX (8.1 and 12.1 /am bands) and Spitzer 
IRAC images exhibit bright nebulous emission coinci- 
dent with the inte rior of the ring-like structure seen in 
the MIPS images ( |Kramer et al. 2003 J0rgensen et al. 
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2006). We first noticed the molecular gas shell when 
we detected high velocity 12 CO in the region and fur- 
ther inspection revealed that the channel maps show a 
clear ring-like morphology coincident with the IR nebu- 
losity, with a velocity structure consistent with that of 
an expanding bubble. The CO emission associated with 
CPS 6 is most clearly seen from Vlsr ~ 7 to 8 km s _1 , 
but emission associated with this shell is seen over a 
wider range of velocities where the ring structure is ob- 
served to slightly diminishes in size as the LSR velocity 
increases (see Figure 16). The 12 CO integrated inten- 



sity (for 6.5 < Vlsr < 8.1 km s ) shows a clear ring 
structure with a radius of about 8' (~ 0.6 pc), slightly 
larger than (and enclosing) the ring seen in the Spitzer 
MIPS 24/im image (see Figure [Tt]) . The p — v diagram 
shows an inverse "U" shape wrtn a minimum LSR ve- 
locity at ~ 5.5 km s _1 and a maximum LSR velocity of 
~ 8.5 km s -1 , where the shell emission blends with the 
cloud emission (Figures 17). This velocity structure is 



consistent with that expected of the backside hemisphere 
of an expanding bubble that has a central velocity off- 
set (towards the blue) with respect to the cloud central 
velocity (similar to the scenario presented in Figure [5]i). 
From the channel maps and p — v diagram we estimate 
V eX p — 3 km s _1 and the central velocity to be at Vlsr 
= 5.8 km s _1 . 

The source IRAS 03382+3145, right at the center of 
CPS 6, is the best candidate driver of the surround- 
ing expanding shell. There is not much information 
on this source, nor there are any estimates of its spec- 
tral type. In the c2d catalog, the point source SSTc2d 
J034122.1+315434, detected in all 2MASS and IRAC 
bands, is the closest source to the position (and lies 
within the position error ellipse) of IRAS 03382+3145, 
making it the most likely counterpart to the IRAS source. 
It is classified as a "star" in the c2d catalog (i.e, with NIR 
SED consistent with a stellar photosphere), and has an 
a = —2.65. Yet, its clear association to the Perseus cloud 
imply that this is not likely a typical main-sequence star 
and it is very probable that this is a young pre-main- 
sequence star. There is no reliable MIR detection in 
the MIPS bands because the so urce is confused wit h the 
bright nebulosity in this region ( |Rebull et aLl|2007| ) 



3.7. CPS 7 

CPS 7 was first noticed when, using our 3D visualiza- 
tion technique used to search for high- velocity features in 
the Perseus cloud, we detected a few positions that show 
blue-shifted CO emission, at velocities significantly dif- 
ferent from the cloud velocity, along the southern edge 
of the Perseus molecular cloud (about 20' southeast of 
the IC 348 cluster). Examining the channel maps we no- 
ticed a wide arc of about 70' in extent, tracing approx- 
imately 3/8 of a circle with a radius of approximately 
30' centered around 03 /l 42 m 18 s , 32°06 / 10 // (J2000), that 
is clearly se para te from the rest of the cloud emission 
(see Figures 18 and 19). We detect faint clumpy emis- 
sion associated with CPS 7 down to Vlsr ~ 5 km s _1 , 
but the arc structure is best seen at Vlsr between 6 
and 7 km s _1 . At Vlsr > 6.8 km s _1 emission as- 
sociated with the IC 348 cluster and CPS 6 overlaps 
with the eastern and western ends of CPS 7, respec- 
tively. Additional emission associated with the L1468 



cloud at Vlsr > 7.5 km s _1 confuses the region, mak- 
ing it impossible to isolate the structure associated with 
CPS 7 at these velocities (Figures 18). The p — v dia- 
gram in Figures [19] shows a clear velocity gradient in the 
CO emission associated with CPS 7 (i.e., the emission at 
angular offset from -22' to -40'). This is consistent with 
the velocity structure expected from part of the backside 
hemisphere of an expanding bubble. From the channel 
maps and p — v diagram we estimate an expansion veloc- 
ity of 3 km s _1 and the central velocity to be at Vlsr 
= 6.0 km s _1 . Although there is no clear IR counter- 
part to the CPS 7 structure observed in CO, it does not 
necessarily mean that it does not exist. This region of 
the cloud is filled with the bright nebulous IR emission 
associated with CPS 5 and CPS 6 to the west and the IC 
348 cluster (and CPS 8) to the east. Thus any faint IR 
emission associated with CPS 7 may be confused with 
the bright emission from these nearby shells. The mass 
region is constrained to the part of the CPS 7 annulus 
that lies west of 03 /l 44 m , in order to avoid emission as- 
sociated with IC 348. 

The most promising candidate driving source for CPS 7 
is IRAS 03390+3158, which lies approximately 2.6' from 
the estimated center of the shell. This source, which by 
chance is coincident with the IR nebulosity associated 
with the CPS 5 dust ring (from HD 278942), is an opti- 
cal variable star known as V633 Per (with an observed 
variation in V from 13.0 to 13.5 mag), and it is listed 
in the 2MASS and c2d point source catalogs. SSTc2d 
J034210. 9+320817, as it is known in the c2d point source 
catalog, is only detected in IRAC bands 1 and 3 as the 
other two IRAC bands did not cover the area. The source 
is also confidently detected at 12/im with IRAS and 24/im 
in the MIPS1 band. Detection of this source at 60/im and 
70/mi with IRAS and Spitzer , respectively, is probably 
hampered by contamination from the surrounding neb- 
ulous emission from the dust ring associated with CPS 
5. The IRAS point source catalog list a detection at 
100/im, but the flux at this wavelength is probably con- 
taminated by the surrounding nebulosity. From the loca- 
tion of this star in a J-H vs. H-K plot, it is apparent that 
its IR colors are consistent with a reddened K5 to M0 



main-sequence star or reddened classical Be star (Lada 
fc Adams|1992| ) . It is highly unlikely that this star would 
be a late type main-sequence star, as its B and V appar- 
ent magnitudes (15.7 and 13.0 to 13.5 mag, respectively) 
and color (5 — V" = 2.2to2.7 mag) imply that it would 
be very close to the Sun (11 to 32 pc) and have an ex- 
tremely unlikely high extinction {Ay = 2.5 to 4.0 mag) 
for a nearby star. It, seems more plausible that this is a 
young late B or early A star with an extinction of about 
6 to 8 mag, residing in or closely to the Perseus cloud at 
a distance of about 250 pc. 

3.8. CPS 8 

CPS 8 is associated with the bright IR nebulosity ob- 



served in the Spitzer images presented by |Muench et al. 
(2007J that lies in the southern part of the IC 348 cluster. 
We detected this shell when we examined the CO chan- 
nel maps and noted that, south of the central part of the 
main IC 348 cluster, the emission at the northern edge of 
the molecular cloud has a curved (arc-like) morphology 
that is significantly brighter than the surrounding CO 
emission and is coincident with the bright border of the 
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IR nebulosity (see Figures 20 and 21 ). The bright curved 
wall of molecular gas emission, wrnch traces an arc that 
is less than half a circle, changes position with LSR veloc- 
ity, consistent with the kinematics of an expanding shell 
(see Figure 20). Molecular gas emission associated with 



this structure is seen from Vlsr ~ 8 km s _i up to ap- 
proximately 10 km s _1 , but the structure is most clearly 
discerned at Vlsr °° 9.3 km s _1 where it shows a peak 
T m 5 of about 30 K (which corresponds to T ex ~ 34 K , 
assuming the CO is optically thick and using Equation 
1). This is consistent with the gas being heated by the 
same external source responsible for the nebulosity. At 
all velocities where CO emission associated with CPS 8 
is detected, but specially at Vlsr < 8.8 km s _1 , there 
is also emission associated with a filament that extends 
northwest of (and it is unrelated to) CPS 8. The p — v 
diagram shows a rough velocity gradient in the arc-like 
CO emission (i.e., increasing Vlsr with increasing dis- 
tance from the center, for angular offsets from 10' to 
20'). Assuming an expanding bubble is responsible for 
this gradient, it is then most likely that the CO arises 
from part of its front hemisphere. Similar to CPS 5 
and 7 (where the CO shell is limited to less than half 
a circle) it is hard to constrain the shell's expansion ve- 
locity and we can only obtain a lower limit. A reason- 
able estimate, using the p — v diagram and fitting the 
expanding bubble model to the channel maps, is an ex- 
pansion velocity of ~ 2.5 km s _1 and a central velocity 
at Vlsr — 9.6 km s _1 . The mass region in CPS 8 is lim- 
ited to the shell area south of 32° ll 7 , in order to avoid 
emission from the filament in the region. 

A circle with a radius of ~ 13' centered on 03 /l 44 m 10 s , 
32°17 / 20 // (J2000) provides a reasonable visual fit to the 
curved wall of high-intensity CO associated with CPS 
8 Within 3' of this circle's center the only known early 
type star capable of illuminating the bright IR nebula is 
the B1III star o Per (at - 1.9' from the center of CPS 
8). It has been typically assumed that this star is not a 
member of the IC 348 cluster nor of the nearby Per OB 
association, mainly due to its peculiar pr oper motions as 
compared to the rest of the association ( Fredrick||1956| 
de Zeeuw et al. 1999) . How ever, the estimated distance 
to this star by Cernis (1993), about 280 ±50 pc, is consis- 
tent with the average distance to t he cluster of 3 00 ±15 pc 



obtained from different methods ( Herbst||2008 ), and it is 
quite possible that o Per is sufficiently "close (within ~ 1 
pc) to the cluster to be responsible for CPS 8. This bright 
source is observed in the optical and the infrared, and it 
is classified as a YSO candidate in the c2d point source 
catalog, where it is named SSTc2d J034419. 1+321718. 

Evidence that o Per might be responsible for CPS 8 
also comes from the fact that extended 8 /im emission 
(in the Spitzer IRAC map of the region) is coincident 
wi th the CO and 24 /im arc-like structures ( |j0rgensen et 



|aO2006j). It is generally thought that the 8 /am emission 
is dominated by PAHs that are heated by far-UV photons 
and o Per is the only source in the region that seems to 
be capable of producing enough UV photons to heat the 
entire arc associated with CPS 8. The candidate drivers 
of CPS 10 and 11 are able to heat their surroundings and 
form small rings that are seen in the IR, but these only 
extend to no more than 2' from their sources' position 
— too small to be responsible for the structure associated 



with CPS 8. Further indirect evidence that o Per is re- 
spons ible for CPS 8 comes from the study of |Sun et al. 
( |2008l ), who derive the far-UV (FUV) field at different 
positions in the IC 348 region by a fitting a PDR model 
to the observed intensity ratios, using lines of CO, 13 CO 
and [CI]. The derived FUV flux at positions coincident 
with CPS 8 is a factor of 3 to 7 higher than the FUV flux 
expected to arise solely from HD 281159 (the source of 
CPS 10). These results imply that an additional source 
of FUV photons is needed to account for the high FUV 
flux in the region, and o Per seems to be the most likely 
source. 

3.9. CPS 9 

This relatively small shell, seen only in the molecu- 
lar line maps, lies approximately 5' north of the eastern 
part of CPS 8 (see Figure |2l]) . It was first noticed when 
we examined in detail the (JO channel maps in the IC 
348 region. The blueshifted (foreground) cap of the ex- 
panding bubble is observed at Vlsr ~ 7.7 — 8.1 km s _1 . 
At larger LSR velocities the CO shows a ring-like struc- 
ture that increases slightly with radius, and reaches a 
diameter of approximately 5' at its largest extent (see 

Figure |22|). We detect CO associated with CPS 9 all 

-l 



the way up to Vlsr ~ 9.9 km s . However, at 
Vlsr > 9 km s _1 only an arc-like structure is observed 
in the southern part of the shell. It appears that CPS 9 
has blown-out of the low-density (northern) edge of the 
cloud. The CO structure associated with CPS 9 has a 
slight elongation in the southwest-northeast direction at 
Vlsr = 8.6 — 9.0 km s _1 . This morphology could be due 
to a lower expansion rate towards the region with higher 
density in the southwestern part of the shell. We need 
higher angular resolution data to confirm this. We find 
no clear IR counterpart to the CPS 9 structure observed 
in CO. However, this does not necessarily mean that it 
does not exist as this region of the cloud is filled with 
the bright nebulous IR emission associated with CPS 8 
and 10, and the CPS 9 IR emission (it it exists) could be 
confused with other IR nebulosity in the region. 

The p — v diagram in Figure [23] shows an asymmet- 
ric "V-like" structure with a minimum LSR velocity at 
~ 7.5 km s _1 and different maxima LSR velocity at the 
two ends: about 9.0 km s _1 at the northwestern (neg- 
ative angular offset) end and approximately 9.5 km s _1 
at the southwestern (positve angular offset) end. This 
asymmetry is due to the fact that there is little or no 
CO emission in the northern edge of CPS 9 at Vlsr > 
9 km s , as indicated above, and possibly also due to 
emission associated to CPS 8 in the southern part of CPS 
9. The velocity structure in the p — v diagram is consis- 
tent with that of the front hemisphere of an expanding 
bubble (similar to that shown in Figure^). From the 
channel maps and p — v diagram we estimate an expan- 
sion velocity of 2 km s _1 and the central velocity to be 
at Vlsr = 9.1 km s -1 . The mass region for CPS 9 is 
constrained not to include regions included in CPS 8 or 
CPS 10. 

A search in SIMBAD for known stars within 65" of 
the center of CPS 9 yields 35 sources. Two possible 
candidates are V* 695 Per and IC 348 LRL 30. The 
former is an M3.75 star t hat is an optical variable with 
detected X-ray emission (Preibisch & Zinnecker 2001) 
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and is classified as a YSO candidate in the c2d cata- 
log (where it is named SSTc2d J034419.2+320735). IC 
348 LRL 30 is a F0 star that has been classified as a 
YSO candidate in the c2d catalog (where it is named 
SSTc2d J034419. 1+320931) and has been also detected 
in X-ray (CXOPZ J034419.1+320931, Preibisch & Zin- 



necker 2001). These two sources are plotted in Figure 23 



where it can be seen that none of these lie at the exact 
center, V* 695 Per is ~ 55" to the south of the center, 
while IC 348 LRL 30 lies - 60" to the north. Out of these 
two sources, IC 348 LRL 30 seems the most likely as it 
the star with the earliest spectral classification. Another 
F type star in the region is IC 348 LRL 28 (at about 
60" from the center), but it appears to be a foreground 
star (jLuhman et al.|2003| ) so we discard it as a candidate 
source. 

3.10. CPS 10 

CPS 10 is associated with the well known source HD 
281159 in IC 348. A circular structure of about 3' in 
diameter around this source is observed in all Spitzer 
IRAC bands and in the MIPS 24/mi and 70 jam images 
flj0rgensen et al ; |2006| [Rebull et al.|2007[ ). The existence 
ol this 1R nebulosity prompted our search for a similar 
structure in the molecular line maps, and we found that 
the CO emission surrounding this source has a ring-like 
morpho logy with part of the eastern region missing (see 



Figures 



24 



and 



25) 



The presumed back (redshifted) cap 
of the bubble is observed at Vlsr ~ 9 km s _1 and the 
size and shape of the CO structure changes slightly with 
decreasing LSR velocity (Figure 24). The p — v diagram 
shows an upside down "V" structure with minimum LSR 
velocity at approximately 7 km s _1 (Figure 25). This ve- 
locity structure is consistent with that of theoack hemi- 
sphere of an expanding bubble. From the channel maps 
and p — v diagram we estimate an expansion velocity of 
2 km s _1 and the central velocity to be at Vlsr = 7.1 
km s _1 . The driving source (HD 281159) is also known 
as SSTc2d J034434.2+320946 in the c2d catalog, where 
it is classified as a YSO candidate. 



3.11. CPS 11 

In the northeastern part of the IC 348 cluster there 
is a bowl-shaped nebulosity with a diameter of about 2' 
clearly detected in all Spitzer IR AC bands and in th e 
MIPS 24/im and 70/im images flMuench et aT]|2007|) . 
Similar to CPS 10, the circular nebulosity observed with 
Spitzer motivated our search for a shell structure in the 
molecular maps in this region. Our CO observations 
barely resolve the structure associated with CPS 11, yet 
a ring-like structure is discernible at LSR velocities near 
7.3 km s _1 , and emission associated with this structure 
is seen from Vlsr ~ 6.7 km s _1 to about 8.5 km s _1 
(see Figure 26). The CO ring structure almost fully sur- 
rounds the circular IR nebulosity (see Figure 27). Similar 
to CPS 10, the p — v diagram of CPS 11 shows an upside 
down "V" structure with minimum and maximum LSR 
velocity at approximately 6.8 and 8.5 km s _1 , respec- 



tively (Figure 27). This velocity structure is consistent 
with that of theoack hemisphere of an expanding bubble. 
From the channel maps and p — v diagram we estimate an 
expansion velocity of 2 km s _1 and the central velocity 
to be at Vlsr — 6.9 km s _1 . 



Inside the nebulosity there are two relatively bright A 
stars, and it is the closest to the center of the nebula (and 
the one that lies ~ 1' from the center of CSP 11), that is 
the most likely driving source; the star named Dust-Bowl 
in SIMBAD. This is an AO star that is also known as IC 
348 LRL 3, 2MASS J03445064+3219067, and SSTc2d 
J034450. 7+321906, and is a known X-ray emitter (C X- 
OPZ J034450.7+321904, [Preibisch fc Zinnecker|[200T] ) . 

3.12. CPS 12 

CPS 12 lies in the northeastern part of the B5 dark 
cloud, at the easternmost edge of the Perseus molecu- 
lar cloud complex. The ring structure, with a radius of 
about 6' is clearly seen in both 12 CO and 13 CO maps 
as well as in the MIPS (160/im) map. The emission in 
the MIPS 24/im and 70/im maps associated with this 
structure is very faint and the ring-like structure is barely 
discernible at these shorter wavelengths. The ring is best 
seen in 12 CO and 13 CO at LSR velocities between 9.5 and 
9.9 km s _1 , but molecular gas emission associated with 
this structure is seen from Vlsr ~ 8 to 10 km s _1 (see 
Figure [28| . The radius of the structure decreases slightly 
as the LSR velocity increases, and it appears that in CPS 
12 we are mostly observing the back side of an expand- 
ing bubble. This is confirmed by the p — v diagram in 
Figure [29] , which shows an upside down "V" structure 
(similar to that shown in Figure [5]i) . Using the channel 
maps and p — v diagram, we estimate an expansion ve- 
locity of 2.5 km s _1 and a central velocity at Vlsr = 8.2 
km s _1 . 

This circular cavity in the mole cular cloud has been re- 
ported before by several authors (|Goldsmith et alT]|1986 
|Langer et al. 1989 Yu et al. 1999p and all agree that 



the star at the center, IRAS 03446+3254, is the source 
responsible for this structure. This source, also known 
as SSTc2d J0 34747. 1+330404, is believed to be a young 

[19861 |Myers et alj | l987 ) 



stellar object (Goldsmith et al. 



and the presence oi a disk is deduced from IR excess and 
Spitzer IRS spectra that show 10 and 20 fim silicate lines 
pres umably from dust in a flatt ened circumstellar struc- 
ture ( jKessler-Silacci et al. 2006). Moreover, the detection 
of [NellJ in the IRS spectrum of this source is believed 
to arise from a disk heated by X -ray radiation from the 
central star (Lahuis et al. 2007[ ). This source, however, 
is not classified as a YSO candidate in the c2d catalog, 
as this area was not covered by the c2d IRAC maps. 

4. ANALYSIS & DISCUSSION 
4.1. Mass, Momentum, Energy of Shells 

To obtain the mass of the shells, we utilized a method 
that uses the 12 CO(1-0) to 13 CO(1-0) ratio to correct 
for the opacity in the 12 CO (1-0) line, as a function of 
velocity. In the general, the line opacity increases the 
closer to the central cloud velocity. Using an optically 
thick line without properly correcting for its velocity- 
dependent opacity will result in an underestimation of 
the mass, momentum, and kinetic energy. 

A brief description of the method is dis cusse d below, 
but for more det ail see Arce et al. ( |2010[ ) and lArce fc 
Goodman (2001). The region used to obtain the mass 
tor each CPS (i.e., the shell mass region) is an annulus 
with the same radius as that of the CPS and a thick- 
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ness slightly wider than the thickness of the shelQ For 
each shell we calculate average spectra of 12 CO (1-0) and 
13 CO (1-0) over the defined shell mass region in order to 
estimate the ratio of 12 CO (1-0) to 13 CO (1-0), hereafter 
denoted #12/13, as a function of velocity. The line ratios 
obtained for each shell mass region were then fitted with 
a second order polynomial. To calculate the shell mass 
at a given position, we directly use the 13 CO emission 
(given by the main beam corrected antenna temperature 
of the line, X^ 3 ) whenever is detected. At large veloci- 
ties away from the central cloud velocity, where the 13 CO 
emission was not reliably detected, we use the 12 CO (1-0) 
emission (X^ 2 ) and the fit to #12/13 (v) to estimate the 
value of X^ 3 at the given velocity and position, using the 
simple equation T^ 3 = T^ b /R 12/ 



Tmb/Riz/isiv)' We then obtain a 



value of the 13 CO o pacity (T13) using equation 1 of Arce 



Goodman (2001), a value of the 13 CO column density 
(-/Vi a ), and subsequently the m ass, using equations 3, and 



4 of Arce & Goodman (2001). We only use spectral data 
that is greater than or equal to twice the rms noise of 
the spectrum. To obtain the total molecular hydrogen 
column density we used the region-dependent values of 
[H2l/[ 13 CQ1 obtaine d for the Perseus molecular cloud by 



Pineda et all (|2QQ8|) and a rat io of 12 CO to 13 CO of 62 
lrom [Langer & Penzias| ( |1993| ). 

We estimate the excitation temperature of the cloud 
gas, T exi for each shell region assuming that the 12 CO(l- 
0) line core is optically thick. We measured the peak 
X^ 2 of each spectrum in the shell, and use Equation 1 to 
obtain a distribution of the ambient cloud T ex values in 
the region of interest. Most shells have T ex from about 
15 to 30 K, and a few exhibit excitation temperatures of 
more than 30 K. We assume a temperature of 25 K for all 
shells, and caution that this may be a lower limit as it is 
very probable (and expected) that the actual gas temper- 
ature is significantly higher than the value of T ex derived 
from the CO (1-0) data. Studies of collimated outflows 
that observe different transitions of the same molecule 
show that temperatures of the entrained molecular gas 
are higher than the excitation temperature of the am- 
bient gas and is not uncommon for the outflo wing gas 
to reach temper atures of 50 to 100 K (e.g., [Hirano fc| 
| Taniguchi||2001[ ). The study by |Beaumont fc" Williams| 
(|2010| ) indicates that in the shells around high mass star 
forming regions they observed with the CO (3-2) line, the 
gas is most likely to have a kinetic temperature of 40 to 
60 K — higher than the excitation temperature obtained 
from the peak CO (3-2) emission. The uncertainties in 
gas temperature, gas opacity, and shell mass region re- 
sult in an uncertainty in the shell mass estimate of about 
a factor of two. 

Using our estimate of the shell mass and expansion 
velocity (see above) we can then obtain the shell mo- 

4 We include an annulus that is slightly wider than the shell's 
thickness as not all emission associated with the shell is included 
within the ring defined by the profile's FWHM (see Figure pi. The 
annulus is widen by increasing (and decreasing) the outer (and in- 
ner) radii by three pixels (i.e., 70") in the CO maps, always avoid- 
ing including emission from unrelated structures. When needed, 
we use different ring radii for different velocities, where the radius 
and thickness for each velocity range is defined with the same pro- 
cedure as described above (using the azimuthally averaged profile 
of the integrated intensity map over that particular velocity range). 



mentum and kinetic energy. These are given by P s heii = 
M shell l/ exp and £ she ii = 0.5M she iiKL 5 respectively, and 
are shown for each shell in Table [2] Figure [30] shows 
a graphic comparison of the relative mass, momentum 
and energy from each CPS. Note, that for most shells 
P s heii and i£ s heii are lower limits as V exp is the minimum 
possible expansion velocity, as discussed above. 

4.2. Shell Shapes and Cloud Geometry 

The distribution of molecular gas associated with the 
shells can be used to constrain the geometry of the ex- 
panding bub ble and to some extent the th ickness of the 
host cloud (Beaumont & Williams 2010). As stated 
above, although all shells in Perseus show evidence of 
expansion, most are only observed over a limited range 
of velocities. For some of the shells only the front 
(blueshifted) hemisphere is seen (CPS 1, 2, and 9), while 
in others only the back (redshifted) hemisphere is seen 
(CPS 6, 10, 11, and 12). Only one source shows blue- 
and red-shifted emission at the center of the circular shell 
(CPS 3), consistent with having front and back caps (al- 
though the p — v diagram o f thi s shell is not fully compat- 
ible with this picture, see § 3.3 for more detail). The four 
largest shells (CPS 4, 5, 7, and 8) show neither front or 
back cap and the CO shells only traces arcs, not complete 
circles. We expect the winds that drive the observed 
shells to be spherical (see below), and we only detect the 
shells in CO where the winds are entraining the cloud; 
we assume that parts of the shells not detected in CO 
(including caps) are regions where there is no molecu- 
lar gas. The sources of these four shells are observed at 
the edge of the cloud and in some cases the sources are 
thought to be in front or behind the cloud. It is therefore 
not surprising that for these large shells the CO does not 
fully trace circles as the bubbles driven by these sources 
are not immersed in the cloud. 

The other eight shells presented in this work are seen 
in projection to be well within the cloud boundaries and 
would be expected to be detectable in CO at all velocities 
if they were fully surrounded by molecular gas material. 
However, as indicated above most of the shells show only 
one hemisphere (or part of it). This indicates that these 
shells lie closer to the cloud's front (shells that only show 
the redshifted or back hemisphere) or back (shells that 
only show the blueshifted or front hemisphere). This 
could happen if either: 1) most shell sources were at 
the front or back edge of a thick (~ 1 pc) cloud; or 2) 
the clouds are thin along the line of sight (comparable 
to the shell radius) and a slight offset in the position 
of the bubble's center with respect to the cloud's center 
results in one hemisphere (the one that is observed in 
CO) entraining most of the gas while the other entrains 
very little or no gas (and is not seen in CO). The former 
scenario seems unlikely as it would require that most 
sources be in the front or back of the cloud. The shell 
sources are thought to be young (~ 1 — 2 Myr, see below) 
and given the low velocity of young stars with respect 
to their cloud (see above), they could not have moved 
much more than 0.3 pc from their birth site. Hence, a 
more likely scenario is that the parent cloud is relatively 
thin along the line of sight (as suggested by Beaumont 
& Williams 2010). 

The medium-size shells (CPS 1, 2, 6, and 12) have a 
radii ranging from about 0.4 to 0.5 pc, and show only 
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one hemisphere. This suggests that the cloud thickness 
(along the line of sight) is not likely to be much more than 
approximately 0.5 pc in the regions close to these shells. 
In the B5 cloud (the host cloud of CPS 12), the major 
and minor axes on the plane of the sky are about 3.5 
and 1.5 pc. From the discussion above, we estimate the 
cloud thickness to be about 15% to 30% of the axes on the 
pl ane of the sky, consistent with the aspect ratios found 



by [Beaumont fc Williams (2010) in a study of bubbles 
in high- mass star forming regions. One would expect 
that many of the small bubbles in Perseus, traced by the 
shells with radii between 0.1 and 0.2 pc, would be fully 
surrounded by molecular gas. Yet, all the small shells in 
IC 348 (CPS 9, 10, and 11) only show the front or the 
back hemisphere. This is probably because these sources 
lie on a relative narrow filament (0.35 pc wide on the 
plane of the sky) that juts out from the northern part 
of the cloud. From the size of the shells, it appears that 
the thickness (along the line of sight) of this filament is 
not much more than 0.2 pc — consistent with the fact 
that we expect the filament's thickness to be similar to 
its width (on the plane of the sky). The only shell that 
shows evidence of possibly being a bubble fully immersed 
in the cloud is CPS 3. This shell has a radius of about 
0.15 pc and lies close to the central part of the cloud 
associated with NGC 1333, where we would expect the 
cloud to be more than 0.3 pc thick. 

4.3. Origin of Shells 

Cavities of different sizes and shapes have been found 
in many clouds forming low- and high-mass stars. How- 
ever, the shells we report here are different from those 
previously reported in the literature. Shells surround- 
ing young or evolved high-mass stars stars are com- 
monly traced by atomic (HI) maps (Chu 2008; Cappa 
et al. 2003, and refereces ther ein), molecular line maps 
Beaumont & Williams 2010) , and observations of the 
ust IK emission jChurchwell et al.||2006"l |2007[ ). Most 
of these have radii of more than 1 pc and it is clear that 
they are produced by stellar winds from one or more 
high-mass stars inside the shell. In nearby regions where 
low- and intermediate-mass stars are forming, cavities 
associated with active and fossil outflows are common. 
One example is NGC 1333, where the approximately 
twenty cavities that have been found are substantially 
smaller (with sizes 0.1 to 0.2 pc), and most have lower 
expansion velocities than the shells reported here ( |Knee| 
' fc Sandel l 2000; Quillen et aLl|2005| . Many of the cav- 
ities in 1NGC l3o3 cannot be associated with a central 
source, and it has been suggested that these are fossil 
outflow lobes — similar to the highly elongated cavities 
found in the Circinus cloud by Bally et al. (1999). On the 
other hand the shells reported here have a circular mor- 
phology and are much larger than typical outflow-driven 
cavities, yet are (in general) smaller than the bubbles as- 
sociated with high-mass stars. Moreover, most (ten out 
of the twelve) of the shell candidate sources in Perseus 
are low and intermediate mass stars (see below), very dif- 
ferent from the well-known bubbles in our galaxy which 
are mostly d riven by high-mass sta rs that produce HII 
regions (e.g., [Deharveng et al.pOlO ). 

The difference between the shells in Perseus and those 
commonly observed near high-mass stars and the out- 
flow cavities that have been reported in other regions of 



active star formation warrants a discussion on the possi- 
ble nature of the COMPLETE Perseus Shells. The fact 
that the majority of the shells coincide with an IR neb- 
ulosity of a similar shape, the existence of a convincing 
candidate driving source (in most cases), and the velocity 
structure of the shells traced by the CO line emission all 
strongly suggest that they are formed by stellar winds 
that drive expanding bubbles in the Perseus molecular 
cloud. The location of the candidate sources, either in- 
side or at the edge of the cloud, as well as their optical-IR 
SED indicate that the sources are likely to be relatively 
young (pre-main sequence) stars. Many of the candidate 
sources show evidence of youth: evidence of a circumstel- 
lar disk, X-ray emission, or optical variability. Except 
in two cases (o Per and HD 278942), all other candidate 
sources with a known spectral type are B5 or later. Even 
for candidate sources without spectral classification we 
can safely deduce from their colors and luminosity that 
they are low or intermediate-mass stars. 

Apart from the Ka emission observed i nside CPS 5 
driven by HD 278942 QRidge et aLl|2QQ6b| ) there is no 
evidence of warm or hot ( > 10 4 K) gas associated with 
any of the other shells reported here. If these expanding 
shells are driven by stellar winds, then we can safely as- 
sume that they are in the momentum conserving (snow- 
plow) phase, as it is typically assumed for winds from 
low-mass star s (see|Norman fc Silk|1980 Levreault|1983 



|McKee| [l989). We can then obtain a rough estimate o: 
the wind mass loss rate ( ra w ) required to produce the 
observed shells using: 



£_shell 



(2) 



where P s heii is the total shell momentum, the wind ve- 
locity is v w , and r w is the wind timescale (the amount 
of time the wind has been active). We obtain the total 
shell momentum from the data (see Table |2| and make 
reasonable assumptions on the values of v w and r w (see 
below) to obtain an estimate of m w . The wind velocity 
is typically assumed to be close to the escape velocity 
of the star, which is about 1 — 4 x 10 2 km s -1 for low - 
and intermediate-mass stars ( Lamers fc Cassinelli|1 999). 
Here we assume v w = 200 km s _1 . The slope of the SED 
between 2 and 20 fim is generally used to estimate the 
ev olution ary stage of a pre-main sequence star (Evans et 
|aT| |2009| ). For most of the candidate sources we obtain 
estimates of a from the c2d point source catalog (see Ta- 
ble |3|. For the most part, the candidate sources appear 
to be in the Class II or Class III stage, which implies tha t 
these have ages of about 1 to 3 Myr (Evans et al. 2009). 
Assuming that the wind has been active tor most of the 
lifetime of the star, then we can assume that r w ~ 1 Myr. 
We use equation 2 to obtain a rough estimate of the wind 
mass loss rate for each source, which ranges from 10 -8 
to 10" 6 M yr _1 (see Table [3). We note that by using 
this formula, and assuming that the winds are spherical, 
we only obtain a lower limit of the wind mass loss rate 
needed to drive the observed shell with P s heii- This is 
because we only detect CO shells where the winds are 
interacting with the molecular gas and for many of the 
shells discussed here the full bubble is not detected in 
CO (see 



4.2). 



We corroborate our estimate of the wind mass loss 
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rates by comparing the measured radius of the shells 
with the expected radius of a wind-driven bubble in the 
momentum-conserving snowplow phase. This is given by: 



3m w v w rl 
2itp 



1/4 



(3) 



where p Q is the density of the interstellar medi um sur 



rounding the bubb le and R\> is the bubble radius ( Lamers 
fc Cassinelli 1999). Using Equation 2, we can then ex- 
press the equation above as: 



i?b = 0.5 pc 



shell 



100 M km s 



1/4 



10 6 yr 



1/4 



5 x 10 4 



stars) or those driven by thermal pressure (i.e., coronal 
winds) could in principle produce circular shells. How- 
ever, these two mechanisms cannot provide a mass loss 
rate that is high enough to explain the observed shells in 
Perseus, as radiation-driven winds from low-mass main 
sequence stars have mass loss rates in the ord er of 10 -12 
to 10 -10 Mq yr -1 ( jLamers fc Cassinelli|l999| ), and coro- 
nal winds produce m ass loss rates of only up to 10 -9 
M yr" 1 flDeCampli|l981 ). 

The inferred properties of the stellar wind responsible 
for the CO shells in Perseus are more consistent with 
tliosf /qf winds detected in pre-main sequence (T Tauri 
cm and \ Herbig Ae/Be) stars. These winds, thought to be 



(4) 

where n Q is the number density of molecular hydrogen 
in the surrounding molecular cloud. Using the value of 
^sheii for each CPS (from Table |5| and the same value of 
the cloud density (n Q = 5 x 10 4 cm -3 ) and wind timescale 
(t w = 1 Myr) in Equation 3, the resulting bubble radius 
(Rb) is within a factor of two for nine of the twelve shells 
in our sample. For the remaining three shells, in one 
of them (CPS 3) the estimated from Equation 4 is 
much larger than the observed CPS radius of 0.14 pc, 
while for two of them (CPS 5 and 7) the estimated R\> is 
much smaller than the observed CPS radius of ~ 2.5 pc. 
CPS 3 is in the central region of the protostellar clus- 
ter NGC 1333, and its driver is probably younger (with 
r w < 10 6 yr) and likely surrounded by higher density 
gas (n Q > 5 x 10 4 cm -3 ) compared to the rest of the 
shell driving sources in our sample (which will result in 
a smaller bubble radius). The two large shells are driven 
by sources that appear to be outside the cloud (or at the 
cloud's edge) where n Q is expected to be significantly less 
than 5 x 10 4 cm -3 (which will result in a larger bubble 
radius) . From our discussion above we can safely assume 
that, for the most part, the shells observed in the Perseus 
molecular cloud are driven by spherical winds with mass 
loss rates between ~ 10 _T — 10 -6 M yr -1 . 

4.3.1. What kind of wind drives the shells? 

Stars drive winds with a range of mass loss rates that 
are produced by different mechanisms. Our results can 
constrain which of these mechanisms are most likely for 
the sources discussed here. In addition to having the 
approximate mass loss rate given above, the winds that 
produced the shells in the Perseus molecular cloud should 
be poorly collimated (i.e., spherical or quasi-spherical) 
winds in order to produce the circular structures ob- 
served in the molecular line maps. Although young em- 
bedded protostars can power winds with mass loss rates 
ranging from 10~ 9 to 10~ 6 yr" 1 ([Lizano et al.||l988[ 

et al. 2006p , consistent with 
winds are typically bi-polar 



|Hartigan et al.||1994| |Bally 

our estimate of ra w " these 



and collima ted (e.g., Arce & Sargent 2005; Stojimirovic 



et al. 2006) and it is highly unlikely that they will pro 



duce the circular morphologies of the CO shells observed 
in our maps. 



(Arce & Sargent 


2006) and young stars 


sequence (Class J 


1) stage may drive out- 



flows with very wide opening angles that could resemble 
spherical winds (see below). The radial winds driven 
by radiation pressure (commonly observed in high-mass 



acc/etion-driven, are known to have mass loss rates of 
up to a few times 10 _T M yr -1 and are believed to 
be spherica l or to have very wide opening ang les p lsmi 



| et al. | [l995l |Edward s et ah] [2006). Recent observations 
by |Edwards et al.| ( |2006l )~of theHe I A10830 profile of a 
sample of accreting T Tauri stars indicate that in many 
sources this line exhibits a P Cygni profile similar to that 
expected from stars with spheri cal stellar winds. Fu rther 
modeling of the line shapes by Kwan et al.l (|2006j) sug- 



gests that about half of the stars power wide-angle polar 
(stellar) winds, and the other half have line shapes that 
are well fitted by a disk wind emerging at a constant an 



gle re lative to the disk surface (see also |Kurosawa et al. 
|2011[ ). These results are consistent with a picture where 
both wide-angle polar winds and disk winds are present 
at the same time in T Tauri st ars, resulting in a quasi- 
spherical wind (see Figure 3 in [Edwards 2009), and the 
wind that will dominate the observed line shape depends 
on the inc lination angle of the system with r espect to the 
observer ( [Kwan et al!]|2006[ |Edwards||2009| . These stel- 
lar and disk winds are accretion-driven winds launched 
through magnetohydrqdynam ic processes in the star-disk 
system ( Ferreira et al.|2006| ). 

We note that although the proposed scenario in which 
the observed Perseus shells are formed by accretion 
driven winds is the most compatible with our results, 
there might be some problems with this picture. First, 
a combination of very wide angle polar winds and disk 
winds does not produce a fully spherical wind. In prin- 
ciple these winds would unlikely be able to produce fully 
circular shells. However, if the direction of the ejection 
axis of the wide angle wind changes on timescales much 
shorter than r w (for example, due to precession of the 
disk, see Terquem et al. 1999), it will allow the wind 
to interact with the cloud over a wider range of angles, 
and could result in a circular CO shell around the source. 
The other problem arises when we try to estimate the re- 
quired reservoir of material in the disk in order to drive 
winds with m w ~ 10 _T — 10 -6 M yr -1 . If these winds 
are active for r w ~ 10 6 yr this means that a total of 0.1 to 
1 M of material will be launched throughout the life- 
time of these winds. Models of accretion-driven winds 
indicate that the ratio of wind mass loss rate (m w ) to 



disk accret ion rate (m acc ) is 0.1 to 0.3 ([Konigl fc Pu-| 
|dritz|2000[ |Shu et al.||2000| |Ferreira et aL|2UUbj ) . Assum- 
ing m w /rn acc ~ 0.3, then disks would need to be more 
massive than what they are typically assumed to be (be- 
tween 0.3 and 3 M ) at the start of the wind phase in 
order to drive winds with high mass loss rates for about 
1 Myr. Recent numerical simulations of embedded sub- 
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solar protostars with disks by Vorobyov (2010) indicate 
that at the end of various runs (representing an age of 
0.7 Myr) the disk to star mass ratio could by as high as 
0.5. If this result scales to the formation of higher mass 
stars, it might then be possible for G, A, and B5 stars to 
have disk masses of about 0.5, 1 and 3 M , respectively, 
at the start of their pre-main sequence phase. 

Our results indicate that out of all the low- and 
intermediate- mass young stars in the Perseus molecular 
cloud only a small fraction drive winds with exception- 
ally high mass loss rates (> 5 x 10 _T M yr _1 ). Out 
of the 12 shells discussed here, only six (CPS 1, 2, 4, 
6, 7, 12) come from presumed low- or intermediate-mass 
stars and require m w > 5 x 10 _T M yr _1 to drive the 
observed shells (see Table |3|. We might be missing a 
few additional shells due to incomplete coverage of the 
cloud or the relative low angular resolution of our CO 
maps. We have no means to estimate how complete is 
our shell sample. Yet, it seems very unlikely that we 
would be missing more than 50% of the existing shells 
produced by winds with fa > 5 x 10 -7 M yr _1 . A 
rough estimate of the total number of young stars in 
Perseus can be obtained using the c2 d catalog of YSO 
candidates for this cloud (Evans et al. 2009). This cat- 
alog shows nearly 400 YSO candidates, but it only lists 
young stars with a detectable infrared excess, and the 
total number of pre-main sequence stars may be slightly 
underestimated. Taking into consideration pre-main se- 
quence stars with no IR excess could increase the number 
of young stars by 10 to 20% (jEvans et al.l|2009[ ). From 
our results (and assuming our shell sample is no more 
than 50% incomplete), we then estimate that roughly 
1 — 3% of the low-mass young stars in Perseus currently 
drive spherical winds with high mass loss rates greater 
than ~ 5 x 10 _T M yr _1 . As stated above pre-main se- 
quence stars power winds with a wide range of mass loss 
rates. Even th ough most have m less than 10 -8 M yr _1 
( |Calvet||l997[ ), winds with up to m ~ 3 x 10~ 7 M yr" 1 
have been detected ( Edwards et al.|[2006 ), and it is not 
implausible that a very small traction of pre-main se- 
quence stars could have winds with mass loss rates as 
high as ~ 10 -6 M yr _1 . Moreover, most of the T Tauri 
stars for which there exists an esti mate of the wind mas s 
loss rate are M or K stars (e.g., |Edwards et al. 2006). 
It is conceivable that slightly more massive (F, A or late 
B) pre-main sequence stars could power winds that have, 
on av erage, higher mass loss rates than those of T Tauri 
stars flNisini et alT||1995 ). 

Observations ol the inner two-thirds of the galactic 
plane by the Spitzer GLIMPSE survey using IRAC im- 
ages revealed approximately 600 ring structures which 
the authors suggest delineate the wall s of bubbles pro- 



duced by w inds from high-mas s stars (Church well et al. 
* ' 2007|) . A recent study by iDeharveng et al.| ( |2010[ ) , 
a sub-sample of the GLIMPSE shells catalogue, 



2006 



using 

showed that about 90% of the shells coincide with ion- 
ized gas from HII regions, traced by radio continuum 
emission. It is possible that the other 10% are formed 
by stars that are too coo l to form detectable HII regions 
(Churchwell et al. 2006| ). In light of the results shown 
here, it could very well be that some of the bubbles 
with no HII regions are powered by late B or A stars 
that are members of the ~ 1 — 3% group of young low- 



and intermediate-mass stars that have exceptionally high 
mass loss rates. 

4.4. Impact of winds on the cloud 

Winds from protostars will inject energy and momen- 
tum into the cloud and may help disrupt their surround - 
ings or drive turbulence in the gas (Nakamura & Li 2007). 
The existence of the shells indicates that stellar winds 
are having a significant impact on the cloud. The most 
obvious effect of the winds is to clear the gas surround- 
ing the young star and pile it up in a ring structure, 
thereby changing the cloud's density distribution. If 
enough gas is accelerated to velocities above the cloud's 
escape speed, the winds could significantly disrupt the 
region. However, from our results it appears that the 
impact from winds to the integrity of the Perseus molec- 
ular cloud complex will be minimal. Using a total mass 
of 7 x 10 3 M for the Perseus cloud complex and an ef- 
fective radius of 7 pc (from the approximate geometrical 
mean of the cloud extent), we calculate the gravitational 
binding energy of the Perseus complex to be ~ 6 x 10 4T 
erg. The total kinetic energy of all expanding shells in 
the Perseus cloud is ~ 7.6 x 10 46 , or approximately 13% 
of the estimated binding energy. Clearly the observed 
shells do not have the energy to potentially unbind the 
entire Perseus molecular cloud complex. However, some 
of the shells are seen at the edge of the cloud and have 
expanding velocities simila r to (or more than ) the escape 



velocity of the region (see Arce et al 



|2010[ ). A fraction 



of the gas can therefore break from the cloud's gravita- 
tional pull, and potentially disrupt the immediate region 
surrounding (within ~ 1 pc) the powering young star. 

The total turbulent energy (E tur b) in the Perseus 
molec ular cloud complex is 1.6 x 10 47 erg ( |Arce et "aL 
|2010[ ). The total kinetic energy of the shells is approx- 
imately one half of E tur ij (see Figure 30), but given the 



uncertainty in the shell mass and expansion velocity es- 
timates the total energy in shells may be comparable to 
the total turbulence energy in Perseus. This suggests 
that stellar winds could potentially help drive the turbu- 
lence in the molecular cloud complex. 

One way to assess whether the winds that produce the 
observed shells can sustain the turbulence in the molec- 
ular cloud is to compare the wind energy injection rate 
into the cloud (E w ) with the turbulent energy dissipation 

rate (L tU rb)' We obtain an approximate value of E w by 
dividing the total kinetic energy of all shells by the wind 
timescale (r w ~ 10 6 yr), and obtain E w ~ 2 x 10 33 erg 
s _1 . If we only consider shells with a confidence level of 
4 or more (see Table EL the wind energy injection rate 
would be about half of the original estimate. Another 
way to estimate t he wind energy injection rate is using 
Equation 3.7 from [McKee] (Il989| : 



En, 



-(M w v w )v 

Y 



(5) 



Here, the total mass loss rate from all the winds in the 
cloud is given by M w . The rms velocity of the turbulent 
motions (i> rms ), is y/S times the velocity dispersion. Us- 
ing the estimate of the wind mass loss rate for each of 
the sources, shown in Table |3j we obtain M w ~ lx 10 -5 
M yr -1 . The average velocity width (FWHM) in 
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Perseus is about 2 km s _1 ( Arce et al.|2010| ), which gives 
u rms ~ 1.5 km s _1 . Assuming v w ~ 2 x 10 2 km s _1 
(see above), we obtain that E w ~ 10 33 erg s _1 , con- 
sistent with our previous estimate. If we only consider 
shells with confidence level 4 or higher, we would obtain 
E w ~ 5 x 10 32 erg s _1 (still consistent, within a factor of 
two, with the previous estimate). The turbulent energy 
dissipation rate is given by L £nr 5 = Eturb/ tdiss, where the 
turbulence dissipation timescale is tdiss — W^ffi which is 
given in terms of the free-fall time (£//). Estimates of the 
value of r] from nu merical simulations of clouds range be- 
tween - 1 and 10 flMcKee|1989l|Mac Low|1999[ ). Assum- 
ing a gas density averaged over the entire cloud complex 
of ~ 10 3 cm -3 and r] = 5, we obtain tdiss ~ 5 x 10 6 yr, 
which results in a turbulent energy dissipation rate of 
approximately 10 33 erg s -1 . The crude estimates of the 
wind energy input rate and turbulent energy dissipation 
rate indicate that these two are approximately the same 
(within a factor of two or less), and our results reveal 
that powerful spherical winds from young stars have the 
potential for driving turbulence in the Perseus molecu- 
lar cloud complex. This conclusion remains even if we 
were to only consider the energy input of shells with a 
confidence level of 4 or more. 

In Perseus, the detected shells are spread throughout 
most of the cloud complex (see Figure |2L yet only a 
small fraction of the cloud (10%, by mass) is currently 
being swept-up in shells. This, however, does not mean 
that they cannot help drive the turbulence in the en- 
tire cloud complex. Clouds are known to be magnetized, 
and it is possible that winds an d shells intera ct with the 
cloud's magnetic field. In fact, Ridge et al.| ( |2006b| ) ar- 
gue that high polarization vectors (from measurements 
of the polarization of background starlight) found to be 
aligned with the warm dust emission from the shell as- 
sociated with HD 278942 (i.e., CPS 5) may be the result 
of the magnetic field in Perseus being swept-up by the 
shell. The interaction between winds and magnetic fields 
may generate large a mplitude Alfven waves that travel 



throu ghout the cloud ( Nakamura & 



2010). These MHD waves could transform 



Li|2007| R 

jistorm tn 
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tions into turbulent motions and help drive turbulence 
in the cloud near and far from the e xpanding shell. 

In their study, |Arce et aL] ( |2010| ) indicate that even 
though outflows have enough power to drive the turbu- 
lence in localized regions of active star formation (with 
sizes of 1 to 4 pc), they lack the necessary energy to feed 
the observed turbulence in the entire Perseus complex. 
The results shown here indicate that, from a pure energy- 
budget perspective, powerful spherical winds from young 
stars can supply the energy needed to maintain turbu- 
lence on a global cloud scale in Perseus. External sources 
of turbulence are probably important during cloud for- 
mation if, for example, clouds are formed through the 
compression of large- scale atomic gas streamers in th e 
interstellar medium (Bal lesteros-Paredes et ah] [1999). 
However, once star formation takes place in the cloud, 
outflows and bubbles can provide enough energy input 
to sustain the cloud turbulence. The radius of the shells 
observed in Perseus range from ~ 0.1 to about 2.8 pc. 
Consequently, shells would be expected to drive turbu- 
lence at a range of scales (similar to outflows) and a dis- 
tinct turbulence driving length would be hard to detect 



see [Arce|[20TT] ). 

5. SUMMARY & CONCLUSIONS 

We detect twelve shells in the COMPLETE molecu- 
lar line maps of the Perseus molecular cloud complex, a 
star forming region typically assumed to harbor low- and 
intermediate-mass stars. The shells, observed in 12 CO(l- 
0) and 13 CO (1-0), have circular or arc-like morphologies 
and range in radius from about 0.1 to 2.8 pc. Most of 
them are coincident with IR nebulosity that have similar 
shape to the CO shells, and for many we can assign a 
candidate source that lies near the center of the circular 
structure. Their velocity-dependent structure is consis- 
tent with the shells being produced by expanding bub- 
bles. All these indicate that the COMPLETE Perseus 
Shells most likely trace the molecular gas that has been 
entrained by spherical (or very wide angle) winds from 
young stars. 

Two out of the twelve shells seem to be powered by 
known high-mass stars, while the rest appear to be driven 
by low- or intermediate-mass pre-main sequence stars. 
One of the high-mass stars is HD 278942, a star that 
is thought to be behind (and possibly associated with) 
the Perseus cloud complex and close enough that its stel- 
lar winds interact with the molecular gas (Ridge et al. 
2006bp . The other is o Per, a star close to 1C 348 that 
has been typically assumed not to be associated with the 
stellar cluster in Perseus. Here we propose that this early 
B star is the most likely source of the relatively powerful 
wind that interacts with the northern part of the Perseus 
cloud and that it lies within 1 pc of IC 348. 

Recent evidence suggests that high-mass stars might be 
impacting the molecular gas in other well-known regions 
that have been typically thought of being predominantly 
low- mass star forming clouds. In Ophiuchus, a ring of 
about 2 pc in radius and approximately centered at the 
position of the B star p Oph is discernible in the dust 
temperature maps, and ap pears to be interac ting with 



the nearby molecular cloud ( Schnee et aLp005| ) . In Tau- 
rus, a cavity of about 3 pc in size is outlined by young 
stars and condensations that look like cometary globules 
( Goldsmith et al.]|2008| ). Although there is not much in- 
formation about the shell- like structure in Taurus and no 
candidate source has been found, it has the makings of 
a shell formed by winds from a high-mass star. It thus 
seem that, contrary to the conventional assumption, low- 
mass star forming regions are not free from the impact 
of winds from high-mass stars. 

Based on corroborating evidence from kinematic signa- 
tures, IR nebulosities, and the identification of candidate 
driving sources, we propose that the shells reported here 
are likely formed by momentum-conserving wind-cloud 
interactions. The observed size and momentum of these 
shells imply stellar mass loss rates of about 10 -8 to 10 -6 
M yr -1 . At least eight require winds with a high mass 
loss rate of m w > 5 x 10 -7 M yr -1 , where two of these 
are the winds the from high- mass stars mentioned above. 
This implies that at least six stars in Perseus drive winds 
that have an unusually high m w (by at least a factor of 
three) compared to the range of mass loss rates measured 
for other pre-main sequence stars. We claim this is not 
unfeasible as no more than 3% of the total young stars 
in Perseus would be required to have such high m w . In 
addition, most of the estimates of wind mass loss rates 
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have been done for low-mass (K and M) T Tauri stars 
and it is plausible that slightly more massive (F, A or 
late B) pre-main sequence stars could power winds that 
have, on average, higher mass loss rates than those of T 
Tauri stars. Although the exact nature of these powerful 
(quasi-spherical) winds is not known, we argue that they 
are likely accretion driven winds. 

We estimate the total kinetic energy of all the shells to 
be approximately one half of the total turbulence energy 
in the Perseus molecular cloud complex, which suggests 
that stellar winds could help drive the turbulence in the 
molecular cloud complex. Furthermore, comparison of 
our estimate of the total energy input rate of all shells 
and the turbulence energy dissipation rate in Perseus in- 
dicate they are similar. Thus, the spherical winds that 
produce the COMPLETE Perseus shells have the poten- 
tial for driving turbulence in the molecular cloud com- 
plex. In Perseus, a comb ination of outflow s in regions 
of active star formation 



6 



et al. 



2010) and shells 



from powerful winds distributed throughout the cloud 
provides more than enough power to maintain the tur- 
bulence throughout the cloud. In principle there is no 
need to invoke large-scale external drivers (e.g., cloud- 
cloud collisions, galactic rotation, etc.) to maintain the 
turbulence in Perseus, at least while the cloud is actively 
forming stars. 

More observations of entire molecular cloud complexes 
using CO line maps (to trace the kinematics of the 
gas) and IR continuum maps (to trace the warm and 
cold dust) with good spatial resolution (~ 0.05 pc or 
less) are needed to study the frequency of shells in 



low-mass star forming regions and their impact on the 
cloud. Further studies should examine whether this is a 
wide-spread phenomenon throughout the Galaxy or the 
Perseus molecular cloud complex is just an exceptional 
case. In addition, large-scale maps of clouds will help 
investigate the frequency at which winds from high-mass 
stars impact nearby low-mass star forming regions. Once 
available, it will be useful to use the molecular line data 



from the JCMT L egacy Survey of the Gould Belt flWard- 
Thompson et al.||2007|) and the IR continuum data from 



the Herschel Gould Belt Survey ( |Andre fc Saraceno|2 005 ) 
to search for shells and study their impact on their parent 
clouds. 
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TABLE 1 
COMPLETE Shells in Perseus 



Shell 


Cloud 


Center 


CO Radius CO Thickness IR Radius 


Candidate 


Name 


Region 


(0:2000, £2000) 


(arcmin/pc) (arcmin/pc) (arcmin/pc) 


Source 



CPS 1 


West of NGC 1333 


03 


27 


29, 


31 


16 


50 


CPS 2 


West of NGC 1333 


03 


27 


36, 


31 


04 


50 


CPS 3 


NGC 1333 


03 


29 


26, 


31 


25 


10 


CPS 4 


East of Bl 


03 


35 


25, 


31 


10 


00 


CPS 5 


Btw. Bl + IC348 


03 


39 


53, 


31 


50 


10 


CPS 6 


L1468 


03 


41 


24, 


31 


54 


10 


CPS 7 


IC348-L1468 


03 


42 


18, 


32 


06 


10 


CPS 8 


IC348 


03 


44 


10, 


32 


17 


20 


CPS 9 


IC348 


03 


44 


19, 


32 


08 


30 


CPS 10 


IC348 


03 


44 


35, 


32 


10 


10 


CPS 11 


IC348 


03 


44 


50, 


32 


18 


10 


CPS 12 


B5 


03 


47 


44, 


33 


03 


50 



6.0 / 0.44 
5.9 / 0.44 
1.9 / 0.14 
16.8 / 1.22 

38.3 / 2.79 
8.0 / 0.58 
31.1 / 2.26 

13.4 / 0.97 
2.5 / 0.18 
2.8 / 0.20 
2.3 / 0.17 
6.0 / 0.44 



3. 


,0 


/ 


0. 


,22 




multiple 


4. 


,7 


/ 


0. 


.34 




multiple 


1 


,5 


/ 


0. 


.11 


1.0 / 0.07 


BD+30 549 


9. 


,4 


/ 


0. 


.68 


11.7 / 0.85 


multiple 


9. 


,1 


/ 


0. 


.66 


36.4 /2.65 


HD 278942 


5, 


,1 


/ 


0. 


.37 


5.5 / 0.40 


IRAS 03382+3145 


7, 


,0 


/ 


0. 


.51 




IRAS 03390+3158 


4, 


,8 


/ 


0. 


.35 


13.0 / 0.95 


omi Per 


1 


,4 


/ 


0. 


.10 




IC 348 LRL 30 


2. 


,5 


/ 


0. 


.18 


1.5 / 0.11 


HD 281159 


1, 


,5 


/ 


0. 


.11 


1.5 / 0.11 


IC 348 LRL 3 


5, 


,0 


/ 


0. 


.36 


5.5 / 0.40 


IRAS03446+3254 



TABLE 2 
Properties of Shells 



Shell 


Mass 


Vexp 


Momentum 


Kinetic Energy 


Name 


[M ] 


[km s- 1 ] 


[M km s _1 ] 


[10 45 erg] 


CPS 1 


48 


2.5 


120 


3.0 


CPS 2 


51 


2.5 


127 


3.2 


CPS 3 


33 


1.2 


37 


0.4 


CPS 4 


66 


5 


330 


16.4 


CPS 5 


53 


6 


315 


18.8 


CPS 6 


124 


3 


371 


11.1 


CPS 7 


92 


3 


277 


8.3 


CPS 8 


151 


2.5 


377 


9.4 


CPS 9 


14 


2 


29 


0.6 


CPS 10 


17 


2 


34 


0.7 


CPS 11 


6 


2 


11 


0.2 


CPS 12 


59 


2.5 


149 


3.7 



TABLE 3 
Properties of Candidate Sources 



Shell 


Candidate 




Position 




Spectral 


c2d 


ex 


rh w a 


Comments 


Name 


Source 




(«2000, <5>2000) 




Type 


Classification 




[10- 7 M yr- 1 ] 




CPS 1 


BD+30 543 b 


03 


27 40.0, 31 15 


40 


F2 


star+dust(MPl) 


-2.47 


6.0 




CPS 2 


multiple 














6.4 


see §[3^2] 


CPS 3 


BD+30 549 


03 


29 19.8, 31 24 


57 


B8 


YSOc_star+dust(MPl) 


-2.33 


1.9 


X-ray source 


CPS 4 


multiple 














15.0 


see pTIl 


CPS 5 


HD 278942 


03 


39 44.7, 31 55 


33 


09.5-B3 


star 


-2.81 


15.8 




CPS 6 


IRAS 03382+3145 


03 


41 22.1, 31 54 


34 




star 


-2.65 


18.6 




CPS 7 


IRAS 03390+3158 


03 


42 10.9, 32 08 


17 




star 


-2.67 


13.9 


variable star 


CPS 8 


Per 


03 


44 19.1, 32 17 


18 


Bl 


YSOc 


-2.66 


18.9 


variable star, X-ray source 


CPS 9 


IC 348 LRL 30 


03 


44 19.1, 32 09 


31 


F0 


YSOc_star+dust(MPl) 


-2.37 


1.5 


X-ray source 


CPS 10 


HD 281159 


03 


44 34.2, 32 09 


46 


B5 


YSOc_star+dust(MPl) 


-2.03 


1.7 




CPS 11 


IC 348 LRL 3 


03 


44 50.7, 32 19 


06 


AO 


star 


-2.65 


0.6 


X-ray source 


CPS 12 


IRAS 03446+3254 


03 


47 47.1, 33 04 


04 




star+dust(MPl) 


-0.14 


7.5 


evidence of disk from IRS data 



a Estimate of minimum wind mass loss rate needed to drive the observed CO shells, using Equation 2, the value of P s heii for each shell given in Table [2] 
r w = 1 Myr, and v w = 200 km s _1 . 

b We assume that the c2d source SSTc2d J032740. 5+311540 is the IR counterpart of BD+30 543. 
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TABLE 4 
CPS Confidence Level Grading 



CPS 


Velocity 


IR 


Circular 


p — V 


Candidate 


Confidence 


Number 


Struct ure? a 


Nebulosity? 15 


Structure? 


Diagram? d 


Source? 63 


Level 


1 


Y 


N 


Y 


Y 


Y 


4 


2 


Y 


N 


Y 


Y 


N 


3 


3 


Y 


Y 


Y 


Y 


Y 


5 


4 


Y 


Y 


Y 


N 


N 


3 


5 


Y 


Y 


Y 


N 


Y 


4 


6 


Y 


Y 


Y 


Y 


Y 


5 


7 


Y 


N 


N 


N 


Y 


2 


8 


Y 


Y 


N 


N 


Y 


3 


9 


Y 


N 


Y 


Y 


Y 


4 


10 


Y 


Y 


Y 


Y 


Y 


5 


11 


Y 


Y 


Y 


Y 


Y 


5 


12 


Y 


Y 


Y 


Y 


Y 


5 



a Does the shell's radius change with velocity as expected for an expanding bubble? 
k Is there IR nebulosity associated with this shell? 

c Does the CO or IR emission associated with this shell show a circular morphology? 
d Is the p — v diagram consistent with that expected for an expanding bubble? 
e Do we find a suitable candidate source? 
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Fig. 1. — Shells in Perseus. The position of each CPS is shown as a dashed circle overlaid on the c2d Spitzer MIPS 1 (24/rni) image (from 
Rebull et al. 2007). 
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Fig. 2. — Shells in the Perseus molecular cloud. Contours show the integrated intensity map of the 13 CO(1-0) emission. Starting contour 
and contour steps are both 3 K km s _1 . These are overlaid on the c2d Spitzer MIPS 1 image. Shells are shown and labeled as in Figure^ 
The names of different clouds in Perseus are labeled in black. 
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Fig. 3. — Three-dimensional volume rendering of the 13 CO (1-0) data around the B5 dark cloud, at the easternmost edge of the Perseus 
molecular cloud complex. Both the color (i.e., rainbow color map) and transparency of each voxel are mapped to the intensity of each pixel 
at each channel. Thus the brightest regions in the channel maps correspond to solid red 3D regions, and the dim emission in the channel 
maps correspond to transparent blue 3D regions. Visible in the upper left portion of the image is CPS 12. This kind of rendering was used 
in |Arce et al.| ( |2010| > to search for high- velocity features. 
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Fig. 4. — Azimuthally averaged radial intensity profiles of three shells. The solid curves show the intensity profiles (from the integrated 
intensity map of each shell). A gaussian fit to the intensity profile was performed on each shell to estimate its radius and thickness. As an 
example, we show the fit to the radial profile of CPS 6 and CPS 7 (dotted line). 
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Fig. 5. — Model of expanding bubble inside a turbulent medium, (a) Integrated intensity map (from -1 to km s 1 ) of model expanding 
bubble. In this model, the bubble center coincides with the center of the cloud that is at a distance of 250 pc, at a Vlsr — km s _1 , and 
it is 1.3 pc thick. The cloud velocity FWHM is 1.2 km s — 1 and the spectral index of the turbulent velocity field ((3) is 2.0. The bubble 
radius, thickness, and expansion velocity are 0.6 pc, 0.25 pc, and 2 km s _1 (see §[3]for a description of the model). This bubble is fully 
immersed in the cloud, (b) p — v diagram of bubble along the direction shown in a. (c) p — v diagram of bubble displaced towards the 
back of the cloud. In this model all parameters are the same as the bubble shown in (a), except the cloud thickness (here is 0.6 pc) and 
the center of the bubble is offset (towards the back) with respect to the cloud's center by 0.4 pc. (d) p — v diagram of bubble displaced 
in position and velocity with respect to the cloud's center. In this model the bubble's center is displaced 0.3 pc towards the front of the 
cloud and the bubble's central velocity is displaced 0.5 km s _1 towards blue velocities with respect to the cloud's central velocity. 




Fig. 6. — Channel maps of 13 CO emission near CPS 1. The number on the upper right corner of each panel indicates the central Vlsr 
of the channel map. Starting contour and contour steps are both 0.5 K. The filled star symbol shows the position of the candidate driving 
source BD+30 543. Dashed circles shows the expected extent, at different radial velocities, of an expanding bubble with a radius, V exp 
and central LSR velocity of 6', 2.5 km s _1 and 8.0 km s -1 , respectively (using the model discussed in §[3J. The location of an unrelated 
filament, CPS 2 and the central part of NGC 1333 are also shown. 
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Fig. 7. — Integrated intensity map and p — v diagram of 13 CO emission near CPS 1. The left panel shows the 13 CO(1-0) integrated 
intensity (for 7.2 < Vlsr < 8-6 km s _1 ) depicting CPS 1. Starting contour and contour steps are both 0.5 K km s -1 . The dashed 
circle shows the approximate extent of CPS 1. The filled star symbol shows the position of the candidate driving source BD+30 543. The 
right panel shows the p — v diagram along the cut shown by the diagonal dashed line in the integrated intensity map. Positions northeast 
(southwest) of the center of CPS 1 are shown as positive (negative) offsets. The approximate extent of the 13 CO emission associated with 
CPS 1 in the p — v diagram is shown as a dashed white rectangle. 
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Fig. 8. — Channel maps of 13 CO emission near CPS 2. The number on the upper right corner of each panel indicates the central Vlsr of 
the channel map. Starting contour and contour steps are 0.6 and 0.5 K, respectively. Dashed circles shows the expected extent, at different 
radial velocities, of an expanding bubble with a radius, V exp and central LSR velocity of 5.9' , 2.5 km s _1 and 8.0 km s _1 , respectively 
(using the model discussed in §[3l. The location of CPS 1 and the central part of NGC 1333 are also shown. 
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Fig. 9. — Integrated intensity map and p — v diagram of 13 CO emission near CPS 2. The left panel shows the 13 CO(1-0) integrated 
intensity (for 6.6 < Vlsr < 8-0 km s _1 ) depicting CPS 2. Starting contour and contour steps are both 0.3 K km s _1 . The dashed circle 
shows the approximate extent of CPS 2. The right panel shows the p — v diagram along the cut shown by the diagonal dashed line in the 
integrated intensity map. Positions northeast (southwest) of the center of CPS 2 are shown as positive (negative) offsets. The approximate 
extent of the 13 CO emission associated with CPS 2 in the p — v diagram is shown as a dashed white rectangle. 
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Fig. 10. — Channel maps of 12 CO emission near CPS 3. The number on the upper right corner of each panel indicates the central Vlsr 
of the channel map. Starting contour and contour steps are 4 and 1 K, respectively. Dashed circles shows the expected extent, at different 
radial velocities, of an expanding bubble with a radius, V exp and central LSR velocity of 1.9', 1.2 km s _1 and 8.1 km s _1 , respectively 
(using the model discussed in §[3l. The filled star symbol shows the position of the candidate driving source BD+30 549. 
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Fig. 11. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 3. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 7.75 < Vlsr < 8-25 km s _1 ) depicting CPS 3, overlaid on the MIPS 24 /im map. Starting contour and contour 
steps are 3.4 and 0.58 K km s _1 , respectively. The dashed circle shows the approximate extent of CPS 3. The filled star symbol shows the 
position of the candidate driving source BD+30 549. The right panel shows the p — v diagram along the cut shown by the diagonal dashed 
line in the integrated intensity map. Positions southeast (northwest) of the center of CPS 3 are shown as positive (negative) offsets. The 
approximate extent of the CO emission associated with CPS 3 in the p — v diagram is shown as a dashed white rectangle. 
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Fig. 12. — Channel maps of 12 CO emission near CPS 4. The number on the upper right corner of each panel indicates the central 
Vlsr of the channel map. Starting contour and contour steps are 0.75 and 1.5 K, respectively. Dashed circles shows the expected extent, 
at different radial velocities, of an expanding bubble with a radius, V e ^ p and central LSR velocity of 16.8', 5 km s _1 and -1.1 km s _1 , 
respectively (using the model discussed in §[3|. The location of the Bl cloud is shown. 
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Fig. 13. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 4. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 1.1 < Vlsr < 2.7 km s _1 ) depicting CPS 4, overlaid on the MIPS 24 /im map. Starting contour and contour 
steps are both 1.5 K km s -1 . The dashed circle shows the approximate extent of the circular IR nebulosity associated with CPS 4. The 
right panel shows the p — v diagram along the cut shown by the diagonal dashed line in the integrated intensity map. Positions northeast 
(southwest) of the center of CPS 4 are shown as positive (negative) offsets. The approximate extent of the CO emission associated with 
CPS 4 in the p — v diagram is shown as a dashed white rectangle. 
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Fig. 14. — Channel maps of 12 CO emission near CPS 5. The number on the upper right corner of each panel indicates the central 
Vlsr of the channel map. Starting contour and contour steps are 0.6 and 2.4 K, respectively. Dashed circles shows the expected extent, 
at different radial velocities, of an expanding bubble with a radius, V exp and central LSR velocity of 38.3', 6 km s — 1 and 6.8 km s — 1 , 
respectively (using the model discussed in § [3}. The filled star symbol shows the position of the candidate driving source HD 278942. 
The location of the main arc associated with CPS 5 is shown in the channel map at Vlsr — 5.5km s _1 . The locations of CO emission 
associated with the IC 348 cluster, the L1468 cloud, CPS 6 and CPS 7 are also indicated. 
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Fig. 15. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 5. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 3.9 < Vlsr < 6.2 km s _1 ) depicting CPS 5, overlaid on the MIPS 24 /im map. Starting contour and contour steps 
are both 2 K km s -1 . The white (black) dashed circle shows the approximate extent of the circular IR nebulosity associated with CPS 5 
(CPS 6). The filled star symbol shows the position of the candidate driving source HD 278942. The right panel shows the p — v diagram 
along the cut shown by the diagonal white dashed line in the integrated intensity map. Positions southwest of the center of CPS 5 are 
shown as negative offsets. The approximate extent of the CO emission associated with CPS 5 in the p — v diagram is shown as a dashed 
white oval. The CO emission associated with the L1468 cloud is also labeled in the p — v diagram. 
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Fig. 16. — Channel maps of 12 CO emission near CPS 6. The number on the upper right corner of each panel indicates the central 
Vlsr of the channel map. Starting contour and contour steps are both 1.5 K. Dashed circles shows the expected extent, at different radial 
velocities, of an expanding bubble with a radius, V^ X p and central LSR velocity of 8 ; , 3 km s — 1 and 5.8 km s — 1 , respectively (using the 
model discussed in §|3}. The filled star symbol shows the position of the candidate driving source IRAS 03382+3145. 
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Fig. 17. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 6. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 6.48 < Vlsr < 8-06 km s _1 ) depicting CPS 6, overlaid on the MIPS 24 /im map. Starting contour and contour 
steps are 1 and 1.5 K km s _1 , respectively. The white dashed circle shows the approximate extent of the CO circular structure associated 
with CPS 6. The filled star symbol shows the position of the candidate driving source IRAS 03382+3145. The right panel shows the p — v 
diagram along the cut shown by the diagonal white dashed line in the integrated intensity map. Positions southeast (northwest) of the 
center of CPS 6 are shown as positive (negative) offsets. The approximate extent of the CO emission associated with CPS 6 in the p — v 
diagram is shown as a dashed white rectangle. 
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Fig. 18. — Channel maps of 12 CO emission near CPS 7. The number on the upper right corner of each panel indicates the central 
Vlsr of the channel map. Starting contour and contour steps are 0.75 and 2 K, respectively. Dashed circles shows the expected extent, 
at different radial velocities, of an expanding bubble with a radius, V exp and central LSR velocity of 31.1/, 3 km s — 1 and 6.0 km s — 1 , 
respectively (using the model discussed in §[3]). The filled star symbol shows the position of the candidate driving source IRAS 03390+3158. 
The location of CO emission associated witn the IC 348 cluster, the L1468 cloud, and CPS 6 is also indicated. 
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Fig. 19. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 7. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 5.52 < Vlsr < 7.11 km s _1 ) depicting CPS 7, overlaid on the MIPS 24 /im map. Starting contour and contour 
steps are both 1 K km s _1 . The white dashed circle shows the size and position of the circle we visually fitted to the CPS 7 CO emission. 
The white dashed oval indicates the position of CPS 6. The filled star symbol shows the position of the candidate driving source IRAS 
03390+3158. The position of the IR nebulosity associated with IC 348 is also labeled. The right panel shows the p — v diagram along the 
cut shown by the thick white dashed line in the integrated intensity map. Positions southwest of the center of CPS 7 are shown as negative 
offsets. The approximate extent of the CO emission associated with CPS 7 (CPS 6) in the p — v diagram is shown as a large (small) dashed 
white rectangle. 
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Fig. 20. — Channel maps of 12 CO emission near CPS 8. The number on the upper right corner of each panel indicates the central Vlsr 
of the channel map. Starting contour and contour steps are 1 and 2 K, respectively. Dashed circles shows the expected extent, at different 
radial velocities, of an expanding bubble with a radius, V exp and central LSR velocity of 13. 4', 2.5 km s _1 and 9.6 km s _1 , respectively 
(using the model discussed in § [3J) . The filled star symbol shows the position of the candidate driving source o Per. 
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Fig. 21. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 8. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 9.27 < Vlsr < 9-65 km s _1 ) depicting CPS 8, overlaid on the MIPS 24 /im map. Starting contour and contour 
steps are 3 and 2 K km s _1 , respectively. The white dashed circle shows the size and position of the circle we visually fitted to the CPS 
8 CO emission. The location of CO emission associated with CPS 9 is shown with a black dashed circle. The filled star symbol shows the 
position of the candidate driving source o Per. The positions of the IR nebulosity associated with CPS 10 and CPS 11 are also shown. The 
right panel shows the p — v diagram along the cut shown by the diagonal thick white dashed line in the integrated intensity map. Positions 
southwest of the center of CPS 8 are shown as negative offsets. The approximate extent of the CO emission associated with CPS 8 in the 
p — v diagram is shown as a dashed white rectangle. 
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Fig. 22. — Channel maps of 12 CO emission near CPS 9. The number on the upper right corner of each panel indicates the central 
Vlsr of the channel map. Starting contour and contour steps are 0.8 and 1.6 K, respectively. Dashed circles shows the expected extent, 
at different radial velocities, of an expanding bubble with a radius, V exp and central LSR velocity of 2.5', 2 km s _1 and 9.1 km s _1 , 
respectively (using the model discussed in § [31 . The filled star symbols in the channel map at Vlsr — 9.8 km s _1 show the position of the 
two candidate sources: V* 695 Per (to the south) and IC 348 LRL 30 (to the north). 
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Fig. 23. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 9. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 8.57 < Vlsr < 9-21 km s _1 ) depicting CPS 9, overlaid on the MIPS 24 /im map. Starting contour and contour 
steps are 2 and 1 K km s _1 , respectively. The white dashed circle shows the extent of the CO emission associated with CPS 9. The location 
of CO emission associated with CPS 9 is shown with a black dashed circle. The filled star symbols show the position of the two candidate 
sources: V* 695 Per (to the south) and IC 348 LRL 30 (to the north). The positions of the IR nebulosity associated with CPS 8, CPS 10 
and CPS 11 are also shown. The right panel shows the p — v diagram along the cut shown by the diagonal thick white dashed line in the 
integrated intensity map. Positions southwest (northeast) of the center of CPS 9 are shown as positive (negative) offsets. The approximate 
extent of the CO emission associated with CPS 9 in the p — v diagram is shown as a dashed (black) rectangle. The CO emission associated 
with CPS 8 is also indicated in the p — v diagram. 
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Fig. 24. — Channel maps of 12 CO emission near CPS 10. The number on the upper right corner of each panel indicates the central 
Vlsr of the channel map. Starting contour and contour steps are 0.8 and 1.6 K, respectively. Dashed circles shows the expected extent, 
at different radial velocities, of an expanding bubble with a radius, V e ^ p and central LSR velocity of 2.6', 2 km s _1 and 7.1 km s _1 , 
respectively (using the model discussed in §|3|. The filled star symbol shows the position of the candidate driving source HD 281159. The 
locations of structures associated with CPS o and CPS 9 are also labeled. 
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Fig. 25. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 10. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 6.8 < Vlsr < 7.8 km s _1 depicting CPS 10, overlaid on the MIPS 24 /im map. Starting contour and contour steps 
are 2.4 and 1 K km s _1 , respectively. The white dashed circle shows the extent of the CO emission associated with CPS 10. The filled star 
symbol shows the position of the candidate driving source HD 281159. The positions of the IR nebulosity associated with CPS 8 and CPS 
11 are also shown. The right panel shows the p — v diagram along the cut shown by the diagonal thick white dashed line in the integrated 
intensity map. Positions southwest (northeast) of the center of CPS 10 are shown as positive (negative) offsets. The approximate extent 
of the CO emission associated with CPS 10 in the p — v diagram is shown as a dashed (white) rectangle. 
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Fig. 26. — Channel maps of 12 CO emission near CPS 11. The number on the upper right corner of each panel indicates the central 
Vlsr of the channel map. Starting contour and contour steps are 0.5 and 1.0 K, respectively. Dashed circles shows the expected extent, 
at different radial velocities, of an expanding bubble with a radius, V exp and central LSR velocity of 2.3', 2 km s _1 and 6.9 km s _1 , 
respectively (using the model discussed in § f3j . The filled star symbol shows the position of the candidate driving source IC 348 LRL 3. 
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Fig. 27. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 11. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 6.7 < Vlsr < 7.4 km s _1 depicting CPS 11, overlaid on the MIPS 24 |Um map. Starting contour and contour 
steps are 0.6 and 0.5 K km s _1 , respectively. The white dashed circle shows the extent of the CO emission associated with CPS 11. The 
filled star symbol shows the position of the candidate driving source IC 348 LRL 3. The positions of the IR nebulosity associated with 
CPS 8 and CPS 10 are also shown. The right panel shows the p — v diagram along the cut shown by the diagonal thick white dashed 
line in the integrated intensity map. Positions northeast (southwest) of the center of CPS 11 are shown as positive (negative) offsets. The 
approximate extent of the CO emission associated with CPS 11 in the p — v diagram is shown as a dashed (white) rectangle. 




Fig. 28. — Channel maps of 12 CO emission near CPS 12. The number on the upper right corner of each panel indicates the central 
Vlsr of the channel map. Starting contour and contour steps are 0.9 and 0.9 K, respectively. Dashed circles shows the expected extent, at 
different radial velocities, of an expanding bubble with a radius, V e ^ p and central LSR velocity of 6 X , 2.5 km s _1 and 8.2 km s _1 , respectively 
(using the model discussed in §[3l. The filled star symbol shows the position of the candidate driving source IRAS 03446+3254. 



COMPLETE Shells in Perseus 



47 




Right Ascension (J2000) Angular Offset 



Fig. 29. — Integrated intensity map and p — v diagram of 12 CO emission near CPS 12. The left panel shows the 12 CO(1-0) integrated 
intensity contours (for 8.4 < Vlsr < 9.7 km s _1 depicting CPS 12, overlaid on the MIPS 24 map. Starting contour and contour 

steps are both 0.8 K km s _1 . The grey dashed circle shows the approximate extent of the CO emission associated with CPS 12. The filled 
star symbol shows the position of the candidate driving source IRAS 03446+3254. The right panel shows the p — v diagram along the cut 
shown by the diagonal thick white dashed line in the integrated intensity map. Positions northwest (southeast) of the center of CPS 12 are 
shown as positive (negative) offsets. The approximate extent of the CO emission associated with CPS 12 in the p — v diagram is shown as 
a dashed (black) rectangle. 
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Fig. 30. — Schematic picture of size, mass, momentum and energy distribution of shells in Perseus. The radius of each circle (and 
position) are proportional to the radius (and location) of the shell in the cloud, while the ring thickness is proportional to the expansion 
velocity (see legend on upper right corner). Shells with a confidence level of 3 or less (from Table [Jl are indicated by a dashed white line. 
Candidate powering sources with a B5 spectral type or later are shown as white stars symbols, while those with earlier spectral type (i.e., 
high-mass stars) are shown as black (filled) star symbols. Candidate sources with no known spectral type (but known a) are shown as red 
stars. The relative mass, momentum and kinetic energy of the shells are shown in the three horizontal bars (where the colors indicate the 
value for each shell). The total outflow mass, momentum and kinetic energy of the molecular outflows in Perseus (from Arce et al. 2010) 
are shown for comparson. 



